AD  -  A  1 32  >30 


COURLIO  ROTOR/AIRFRAMC  VIRRATION  ANAlVSIS(u)  UNITED 
TCCHNOLOOIIS  CORA  STRATFORD  CT  SIKORSKY  AIRCRAFT  OIV 
R  SORHCR  IT  AL .  NOV  >3  NASA-CR-SMS  NASI-  (ION 

_ F/0  1/3 


V 3 


unclassified 


f 


Itt 

|SO 

■EgSj 

t  m 

rm 

lT  «3£ 

fc,  l,IS 

|ii| 

t  li£  1 

*>I.L 

wSi 

I 

L£ 

mm 


MICROCOPY  RESOLUTION  TEST  chart 

national  BuSfAu  or  STANDARDS-  >963  -A 


J 


1  Report  No. 

NASA  CR-3582 


4  Title  and  Subtitle 


COUPLED  ROTOR/AIRFRAME  VIBRATION  ANALYSIS 


3.  Recipient  s  Catalog  No 


5  Report  Date 

November  l 98 2 


6.  Performing  Organization  Code 


7  Author(s) 

Robert  Sopher,  R.E,  Studwell,  S.  Cassarino,  and 
S,B  R.  Kottapalli 


9  Performing  Organization  Name  and  Address 

United  Technologies  Corporation 
Sikorsky  Aircraft  Division 
North  Main  Street 

Stratford,  CT  06602 _ 

12  Sponsoring  Agency  Name  and  Address 

National  Aeronautics  and  Space  Administration, 
Washington,  D.C.  20546 


8  Performing  Organization  Report  No 


10  Work  Unit  No 


1 1  Contract  or  Grant  No 

NASl-16058 

13  Type  of  Report  and  Period  Covered 

Contractor  Report 

14  Sponsoring  Agency  Code 


is.  Supplementary  Notes  Langley  Technical  Monitor:  R.L,  Tomaine 

The  contract  research  effort  which  has  led  to  the  results  of  this  report  was  finan¬ 
cially  supported  by  the  Structures  Laboratory,  AVRADCOM  Research  and  Technology 
Laboratories.  Final  Report. 


16  Abstract 

This  report  describes  a  Coupled  Rotor/Airframe  Vibration  Analysis  which  was  developed 
as  a  design  tool  for  predicting  helicopter  vibrations,  and  a  research  tool  to  quantify 
the  effects  of  structural  properties,  aerodynamic  interactions  and  vibration  reduc¬ 
tion  devices  on  vehicle  vibration  levels.  The  analysis  consists  of  a  Base  Program 
utilizing  an  impedance  matching  technique  to  represent  the  coupled  rotor/airframe 
dynamics  of  the  system  supported  by  inputs  from  several  External  Programs  supplying 
sophisticated  rotor  and  airframe  aerodynamic  and  structural  dynamic  representation. 

The  theoretical  background,  computer  program  capabilities  and  limited  correlation 
results  are  presented  in  this  report.  Correction  results  using  scale  model  wind 
tunnel  results  show  that  the  analysis  can  adequately  predict  trends  of  vibration 
variations  with  airspeed  and  higher  harmonic  control  effects.  Predictions  of  absol- 
utevalues  of  vibration  levels  were  found  to  be  very  sensitive  to  modal  character¬ 
istics  and  results  were  not  representative  of  measured  values. 


17  Key  Words  (Suggested  by  Author(sl) 

Coupled  rotor/airframe  Impedance 


Helicopters 
Vibrations 
Rotor  dynamics 
Airframe  dynamics 


19  Security  Clastif  (of  this  report! 

Unclassified 


matching 


18.  Distribution  Statement 

Unclassified  -  Unlimited 


Subject  Category  05 


20  Security  Classif  (of  this  page) 

Unclassified 


21  No  of  P»9M 


For  sale  by  the  National  Technical  Information  Service  Springfield  Virginia  22161 


NA5M.dnQley.  1982 


NASA  Contractor  Report  3582 


Coupled  Rotor/ Airframe 
Vibration  Analysis 


Robert  Sopher,  R.  E.  Studwell, 

S.  Cassarino,  and  S.  B.  R.  Kottapalli 

United  Technologies  Corporation 
Stratford,  Connecticut 


Prepared  for  Structures  Laboratory 
AVRADCOM  Research  and  Technology  Laboratories 
and  NASA  Langley  Research  Center 
under  Contract  NAS  1-1 6058 


NASA 

National  Aeronautics 
and  Space  Administration 

Scientific  and  Technical 
Information  Branch 


FOREWORD 


This  report  describes  work  performed  by  the  Sikorsky  Division  of 
the  United  Technologies  Corporation  under  Contract  NAS1-16058  for 
the  Structures  Laboratory  of  the  U.S.  Army  Research  and  Technology 
Laboratories  (AVRADCOM). 


The  Technical  Monitor  was  Mr.  R.L.  Tomaine. 


TABLE  OF  CONTENTS 


PAGE  NO. 

Forward .  iii 

List  of  Illustrations .  v 

List  of  Tables .  ix 

Summary .  1 

Introduction .  2 

List  of  Symbols .  5 

Overview  of  the  Coupled  Rotor/Airframe 

Vibration  Analysis .  17 

Substructure  Assembly  Method .  20 

Solution  Modes .  32 

Assembly  of  Rotor  Impedance  Matrix .  42 

Base  Program  Substructures .  48 

Higher  Harmonic  Control  for  Vibration  Reduction .  67 

Sample  Applications .  76 

External  Programs  in  the  Program  Package .  92 

External  Program  Applications.... .  114 

User  Documentation .  120 

Correlation  of  Theory  and  Test .  123 

Concluding  Remarks .  158 

Appendixes 

A  -  Illustration  of  System  Synthesis .  160 

B  -  Substructure  Equations .  169 


References 


193 


LIST  OF  ILLUSTRATIONS 


FIGURE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 


NUMBER 


PAGE  NO. 


Coupled  Rotor/Airframe  Vibration  Analysis 

Data  Flows .  19 

Displacements  and  Forces  Acting  on  a 

Connection  Node  of  a  Substructure .  22 

Local  and  Global  Axes .  26 

Elements  of  Rotation  Transformation 

Matrix  [E] .  30 

Fixed  Absorber  FA1 .  50 

Generalized  Force  GF1 .  52 

Uniform  Beam  Segment  BM1 .  53 

Horizontal  Linear  Bifilar  BF1 .  58 

Vertical  Linear  Bifilar  BF2 .  61 

Basic  Nodal  Isolator .  63 

Vibration  Isolator  IS1 .  65 

Schematic  Showing  Nodal  Isolator  with 

Second  Antiresonant  Bar  Weight .  66 

Bifilar  Mass  Amplitude  Variation  with 

Airspeed .  81 

Normal  Modes  of  a  Ten  Element  Uniform 

Beam  with  Free  Ends .  83 

Variation  with  Time  of  the  Displacement  of 

a  Two-mass  System  to  which  an  Oscillating 

Force  is  Applied .  85 

Frequency  Response  of  a  Simple  Nodal 

Isolator .  86 

Configurations  Used  to  Test  I SI  Com¬ 
ponent .  87 


v 


1 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


FIGURE 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


NUMBER  PAGE  NO. 

Six  Degree-of-Freedom  Isolator .  89 

Frequency  Response  of  Six  Degree-of- 
Freedom  Isolator  at  a  Fuselage  Attach¬ 
ment .  90 

Six  Degree-of-Freedom  Isolator  Vertical 

Responses  at  Fuselage  Attachment  Points .  91 

Effects  on  Rotor  Hub  Loads  of  Variations 

in  Third  Fixed  System  Mode  Frequency  .  93 

Effects  on  Rotor  Hub  Accelerations  of 
Variations  in  Third  Fixed  System  Mode 
Frequency  .  94 

Computer  Wake  Representation-Undistorted 

Wake .  98 

Approximations  to  Lift  Curves  in  F389 

Program .  100 

Calculation  Sequences  for  the  Coupled 

G400/F389  Programs .  101 

General  Arrangement  of  Lifting  Section 

Vortex  Lattice  Network .  104 

APT/Plot  Program  Fuselage  View  from  Below. . .  107 

APT/Plot  Program  Fuselage  View  from  Above...  107 

APT/Plot  Program  View  Enlargement .  108 

Empennage  and  Wake  Lattice  Geometry  in 

RIEVA  Analysis .  110 

An  Example  of  Splitting  of  Rotor  Wake  into 

Two  Parts .  112 

Stabilizer  Airloads  Time  History  in  the 

RIEVA  Analysis .  118 

ARES  Model  in  NASA  Langley  Transonic 

Dynamics  Tunnel  with  Model  Test  Blades .  125 

Sensitivity  of  Rotor  Impedance  to  Hub 
Perturbation  Amplitude .  130 


vi 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


FIGURE  NUMBER  PAGE  NO. 

35  Sensitivity  of  to  HHC  Perturbation 

Amplitude .  130 

36  Correlation  of  Fuselage  Accelerations .  132 

37  Variations  of  Blade  Moments  with  Advance 

Ratio .  134 

38  Radial  Distributions  of  Blade  Moments 

at  0.3  Advance  Ratio .  13  5 

39  Correlations  of  Blade  Moments  at  0.3 
Advance  Ratio  with  Test  Data  from 

Reference  20 .  136 

40  Blade  Moment  Harmonic  Correlation .  137 

41  Blade  Moment  Harmonic  Correlation  with 

ACR  Test .  138 

42  Comparison  of  Variations  with  Advance 

Ratio  of  Predicted  Hub  Loads  and  Measured 
Balance  Loads .  140 

43  Effects  of  Hub  Displacement  on  Blade 

Loads .  142 

44  Effects  of  Hub  Displacement  on  Hub 

Loads .  143 

45  Effects  of  Specified  HHC  on  Amplitudes  of 

Airframe  Vibrations  -  Three  Mode  Solution. . .  145 

46  Effects  of  Specified  HHC  on  Amplitudes  of 

Airframe  Vibrations  -  Nine  Mode  Solution. . . .  146 

47  Effects  of  Specified  HHC  on  Longitudinal 
Vibration  Vectors  -  Three  Mode  Solutions....  147 

48  Effects  of  Specified  HHC  on  Measured 

Blade  Moment  Distributions .  148 

49  Effects  of  Specified  HHC  on  Harmonics  of 

Measured  Blade  Moments .  149 

50  Effects  of  Specified  HHC  on  Predicted 

Blade  Moment  Distributions .  150 

vii 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


FIGURE  NUMBER  PAGE  NO. 

51  Effects  of  Specified  HHC  on  Hub  Forces .  152 

52  Effects  of  Optimal  HHC  on  Fuselage 

Vibrations .  153 

53  Effects  of  Optimal  HHC  on  Blade  Moments .  154 

54  Effects  of  Optimal  HHC  on  Variations  with 
Advance  Ratio  of  Blade  Moments  at  One 

Blade  Station .  155 

55  Effects  of  Optimal  HHC  on  Harmonics  of 
Blade  Moments  at  One  Advance  Ratio  and 

Blade  Station .  156 

56  Effects  of  Optimal  HHC  on  Comparisons  of 
Predicted  Hub  Loads  and  Measured  Balance 

Loads .  157 

B1  Matrices  for  Fixed  Absorber  FA1 .  170 

B2  Mass  Matrix  of  Uniform  Beam  Segment .  173 

B3  Stiffness  Matrix  of  Uniform  Beam  Segment. . . .  174 

B4  Non-Zero  Submatrices  for  Horizontal 

Linear  Bifilar  BF1 .  179 

B5  Mass  Matrix  for  Vertical  Linear  Bifilar 

BF2 .  182 

B6  Damping  Matrix  for  Vertical  Linear 

Bifilar .  183 

B7  Stiffness  Matrix  for  Vertical  Linear 

Bifilar .  184 

B8  Non-Zero  Mass  Submatrices  for  Isolator 

IS1 .  187 

B9  Non-Zero  Stiffness  or  Damping  Submatrices 

for  Isolator  IS1 .  188 

viii 


LIST  OF  TABLES 


PAGE 


TABLE 

NUMBER  NO . 

1  ATTRIBUTES  OF  TARGET  VIBRATION  ANALYSIS .  3 

2  BASE  PROGRAM  SOLUTION  MODES  AND  OPERA¬ 
TIONAL  MODES .  33 

3  BASE  PROGRAM  SUBSTRUCTURES .  49 

4  AEROELASTIC  ROTOR  RE2  COORDINATE  LABELS .  55 

5  CAPABILITIES  OF  G400  ROTOR  AEROELASTIC 

ANALYSIS .  68 

6  BASE  PROGRAM  TEST  CASES  SUMMARY .  77 

7  INPUT  CHECK  CASES  FOR  EXTERNAL  PROGRAMS .  115 

8  ARES  AIRFRAME  MODAL  ATTRIBUTES .  126 

9  CALCULATED  BLADE  NATURAL  FREQUENCIES  FOR 

THE  MODEL  ROTOR  (PER  REV) .  129 

10  COMPARISON  OF  ARES  AIRFRAME  AND  MODEL 

ROTOR  BLADE  NATURAL  FREQUENCIES .  129 


ix 


SUMMARY 


This  report  describes  a  Coupled  Rotor/Airframe  Vibration  Analysis, 
which  was  developed  as  a  design  tool  for  predicting  helicopter 
vibrations,  and  a  research  tool  to  quantify  the  effects  of  struc¬ 
tural  properties,  aerodynamic  interactions  and  vibration  reduction 
devices  on  vehicle  vibration  levels.  The  analysis  consists  of  a 
Base  Program  utilizing  an  impedance  matching  technigue  to  repre¬ 
sent  the  coupled  rotor/airframe  dynamics  of  the  system  supported 
by  inputs  from  several  External  Programs  supplying  sophisticated 
rotor  and  airframe  aerodynamic  and  structural  dynamic  representa¬ 
tion.  The  Base  Program  has  a  stand-alone  capability  to  provide 
assessment  of  the  effects  of  changes  in  airframe  dynamics  and 
vibration  reduction  devices  on  vehicle  vibration  levels.  The  Base 
Program  is  not  restricted  in  input  support  from  any  particular 
programs  and  can  be  executed  interactively  for  many  pertinent 
applications..  The  theoretical  background,  computer  program 
capabilities  and  limited  correlation  results  are  presented  in  this 
report.  Correlation  results  using  scale  model  wind  tunnel  results 
show  that  the  analysis  can  adequately  predict  trends  of  vibration 
variations  with  airspeed  and  higher  harmonic  control  effects. 
Predictions  of  absolute  values  of  vibration  levels  were  found  to 
be  very  sensitive  to  modal  characteristics  and  results  were  not 
representative  of  measured  values. 
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INTRODUCTION 


The  reduction  of  helicopter  vibrations  has  become  increasingly 
important  in  recent  years.  Motivations  are  improvements  in 
reliability  and  maintainability,  passenger  and  crew  comfort, 
visual  acuity  and  increased  equipment  protection.  The  need  to 
reduce  vibration  levels  has  influenced  the  improvement  of  vi¬ 
bration  analysis  methodology.  In  order  to  provide  an  analyses 
capable  of  accurate  vibration  prediction,  the  physical,  mathe¬ 
matical  and  software  attributes  listed  in  Table  1  must  be  in¬ 
cluded.  A  review  of  the  literature  (see  References  1  to  3 )  showed 
that  published  analyses  do  not  possess  a  complete  set  of  these 
desired  attributes. 

From  a  physical  standpoint,  the  adequate  representation  of  the 
vibration  model  requires  an  adequate  simulation  of  the  sources  of 
vibration,  an  elastic  fuselage,  and  the  elements  used  to  control 
the  vibration.  Sources  of  vibration  are  the  air  forces  acting  on 
the  rotor,  fuselage,  and  fixed  surfaces  of  the  helicopter. 
Vibration  control  methods  include  reduction  of  vibratory  loads  at 
the  sources,  modification  of  response  dynamics,  and  isolation  of 
vibration.  Limitations  occur  in  one  or  more  of  these  categories 
and  a  complete  physical  model  does  not  appear  to  be  available  with 
the  scope  desired  for  design  studies. 

Regarding  the  mathematical  methodology,  the  majority  of  methods 
are  limited  in  their  flexibility  to  represent  configurations  and 
in  their  expandability  to  new  configurations  since  they  are  based 
on  hand-derived  equations  for  the  coupled  dynamical  system,  in 
contrast  to  automated  coupling. 

From  a  software  standpoint,  there  are  deficiencies  in  organiza¬ 
tion,  and  internal  documentation.  Comprehensive  modeling  capa¬ 
bility  and  efficiency  require  good  software  to  facilitate  the 
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TABLE  1 


ATTRIBUTES  OF  TARGET  VIBRATION  ANALYSIS 


Physical 

.  Comprehensive  rotor  model 

Articulated,  hingeless,  teetering 
Fixed  system 
Modal 
Stick 
Hybrid 
Drive  Train 
Test  Data 

Vibration  Control  Devices 

Absorbers  (hub,  blade,  airframe) 

.  Isolators 
Higher  harmonic  control 
Other  vibration  sources 
Tail  rotor 
Empennage 
Aerodynamic 

.  Rotor  induced  flow  field 
Stall 

Flight  modes 

Hover,  forward  flight,  transition 


Mathematical 

General  Configuration  modeling  ability 

Steady  state  forced  response,  eigensolution,  time  history 
Parameter  variation 


Software 

Efficiency 

User-oriented  features  -  input,  plotting 
Clear  Architecture 
Code  clarity 
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addition  of  new  components ,  and  the  latitude  to  exercise  and 
improve  solution  modes  and  algorithms.  The  majority  of  programs 
reviewed  were  designed  before  it  was  recognized  that  the  accom¬ 
plishment  of  good  software  is  important  for  flexibility  and  ease 
of  growth. 

In  order  to  provide  the  desired  attributes  a  substructure  assembly 
method  was  utilized  to  simultaneously  provide  an  efficient 
computer  code  while  retaining  the  required  high  level  of  aero¬ 
dynamic  and  dynamical  representations.  The  substructure  assembly 
method  incorporates  a  stand-alone  Base  Program  representing  the 
coupled  rotor/airframe  system  supported  by  External  Programs  which 
provide  the  necessary  aerodynamic  sophistication. 

This  report  describes  the  basis  of  the  coupled  rotor/airframe 
analysis,  the  capabilities  and  limitations  of  the  analysis,  and 
theoretical  and  practical  applications. 
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LIST  OF  SYMBOLS 


Matrix  in  complex  eigensolution  equation  (see 
equation  (36)),  2N  x  2N 

Inverse  of  system  impedance  matrix  (see  equation 
(86)  ) ,  2N  x  2N 

Matrix  in  complex  eigensolution  equation  (see 
equation  (37)),  2N  x  2N 

Matrix  in  Newmark-Beta  method  (see  equation 
(53)),  N  x  N 

Cosine 

Independent  coordinate  system  damping  matrix, 

N  x  N 

Dependent  coordinate  system  damping  matrix,  NQ  x 
ND 

Damping  matrix  for  i-th  substructure 

Matrix  in  Newmark-Beta  method  (see  equation 
(48) ),  N  x  N 

Rotation  matrix  of  direction  cosines  (see 
Figure  4),  3x3 


Vector  of  external  forces  for  i-th  sub¬ 
structure,  representing  forces  and  moments. 


Hll'  H12 ' 
H21 '  H22 


Components  of  force  applied  to  a  substructure 
connection  node  in  x.^ ,  y  ,  and  directions, 
respectively  (see  figure  (2)),  N  (lb). 


Independent  coordinate  system  force  vector, 
N  x  1 

Matrix  in  Newmark-Beta  method  (see  equation 
(54)),  N  x  N 


Dependent  coordinate  system  force  vector,  ND  x  1 


Independent  coordinate  system  force  vector  at 

time  t  in  Newmark-Beta  method,  N  x  1 
n 

Acceleration  of  gravity,  m/s2  (ft/s2) 

Time  interval  in  Newmark-Beta  method  (see  equa¬ 
tion  (43)),  s 


Matrix  relating  rotor  hub  reaction  forces  to  HHC 
control  perturbations  (see  equation  (80)),  12  x  6 


hllH'  h12H ' 
h21H'  h22H 


Submatrices  of  matrix  H  (see  equation  (81)), 
N  x  3 


Independent  coordinate  system  matrix  introducing 
effects  of  HHC  control  perturbations  on  system 
vibrations  (see  equations  (84)  and  (87)),  2N  x  6 


Submatrices  of  matrix  H  (see  equation  (89)), 
N  x  3 


H11D'  H12D'  Dependent  coordinate  system  submatrices  mtro- 

H21D  H22D  ducing  effects  of  HHC  control  perturbations  on 

dependent  coordinate  system  reaction  loads  (see 
equations  (82)  and  (83)),  x  3 

I  Identity  matrix 

J  Performance  index  for  vibration  minimization  (see 

equation  ( 90 ) ) 

K  Independent  coordinate  system  stiffness  matrix,  N 

x  N 


“d 

[k](l) 

M 


Dependent  coordinate  system  stiffness  matrix, 
x  Nd 

Stiffness  matrix  for  i-th  substructure 

Independent  coordinate  system  mass  matrix, 

N  x  N 


"d 


m 

m 

m 


xl 

yi 

zl 


(ml 


(i  ) 


Number  of  normal  modes  in  a  substructure  modal 
representation 

Dependent  coordinate  system  mass  matrix,  x 

Components  of  moment  applied  to  a  substructure 
connection  node  in  x^f  y^,  and  z ^  directions 
respectively,  N-m  (lb-ft) 

Mass  matrix  for  the  i-the  substructure 


n 


integer  step  numbers  in  time-history  solution 
(see  equation  (43)) 
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Index  denoting  n-th  bifilar  mass 


N 


N 


D 


P 


<3 


Q 


{r} 


(i) 


rmxl 

rmyl 

rmzl 


Total  number  of  independent  coordinates  in  the 
assembled  system 

Number  of  bifilar  arms 

Number  or  rotor  blades 

Total  number  of  coordinates  in  the  dependent 
coordinate  vector  (obtained  by  summing  the 
coordinates  for  all  substructures) 

Matrix  in  Newmark-Beta  method  (see  equation 
(49  )  ) ,  N  x  N 

Modal  coordinate 

Matrix  in  Newmark-Beta  method  (see  equation 
(50)  ),  N  x  N 

Vector  of  reaction  forces  for  i-th  sub¬ 
structure,  representing  force  and  moment 
reactions 

Components  of  reaction  force  applied  to  a  sub 
substructure  connection  node  in  x^,  y^ ,  and  z^ 
directions,  respectively  (see  Figure  (2)), 

N  (lb) . 

Components  of  reaction  moment  applied  to  a 
substructure  connection  node  in  x^,  y^ ,  and  z^ 
directions,  respectively  (see  Figure  (2)), 

N-m  (lb-ft) 
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R 


Matrix  in  Newmark-Beta  method  (see  equation 
(51)  ),  N  x  N 


S  Scalar  coefficient  of  n-th  harmonic  of  stress  or 

n 

pushrod  load  on  a  rotor  blade,  N/m2  ( lb/in2 ) 

(S  )  Scalar  coefficient  of  n-th  harmonic  of  stress  or 

no 

push  rod  load  on  a  rotor  blade  for  no  hub  motion 
and  no  HHC  control  input  (see  equation  (94)), 
N/m2  ( ib/in2 ) 

[S^]  Matrix  introducing  the  effects  of  hub  displace¬ 

ment  on  n-th  harmonic  coefficient  of  blade 
stress  or  push  rod  load  (see  equation  (94)), 

1  x  12 

[S  ]  Matrix  introducing  the  effects  of  HHC  perturba- 

n0 

tions  on  n-th  harmonic  coefficient  of  stress  or 
push  rod  load  (see  equation  (94)) 

t  time  at  n-th  integration  step  (see  equation 

(43)),  s 

T  Transfer  matrix  relating  HHC  inputs  to  vibratory 

response  of  the  system,  (see  equations  (84)  and 
(87  )  ) ,  2N  x  6 

u  Xj^  displacement  of  a  connection  node  (see  Figure 

(2)),  m  (ft) 
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(j  ,  Up  Components  of  airflow  relative  to  a  point  on  a 

rotor  blade,  parallel  and  perpendicular  to  the 
rotor  rotation  plane,  respectively,  m/s  (ft/s) 

v  Yp  displacement  of  a  connection  node  (see  Figure 

2),  m  (ft) 

w  Zp  displacement  of  a  connection  node  (see  Figure 

2),  m  (ft) 

Wp  Component  of  airflow  normal  to  the  median  surface 

of  the  rotor  blade  excluding  the  rotor  induced 
flow  (see  equation  (102)),  m/s  (ft/s) 

Response  vector  weighting  matrix  in  expression 
for  vibration  reduction  performance  index  (see 
equation  (90)),  x  NQ 

WQ  HHC  control  angle  weighting  matrix  in  scalar 

performance  index  expression  (see  equation  (90)), 
3x3 

{x} ( 1  *  Vector  of  coordinates  for  i-th  substructure 

X  Vector  of  independent  coordinates  for  the  as¬ 

sembled  system  (=Xp),  N  x  1 

X  Coefficient  vector  of  cosine  harmonic  term  in 

c 

harmonic  expansion  of  X  vector  (see  equation 
(24)),  N  x  1 

X  Coefficient  vectoi'  of  cosine  harmonic  term  in 

co 

i armonic  expansion  of  Xj  vector  for  no  HHC 
control  input  (see  equations  (85)  and  (86)), 

N  x  1 
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D 


<Vc 


<Vs 


II 
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(XT  ) 

'  I  'n 


Vector  of  coordinates  for  the  dependent  coord¬ 
inate  system,  x  1 

Coefficient  vector  of  cosine  harmonic  term  in 
harmonic  expansion  of  vector,  ND  x  1 

Coefficient  vector  of  sine  harmonic  term  in 
harmonic  expansion  of  vector,  ND  x  1 

Vector  of  independent  coordinates  for  the  system 
(=  X) ,  N  x  1 

Vector  of  independent  coordinates  for  system 
after  transformation  (see  equation  (14)) 

Vector  of  independent  coordinates  for  the  system 
after  transformation  P2  (=Xj=X)  (see  equation 
(15)),  N  x  1 

Vector  of  independent  coordinates  at  n-th  time 
step  in  time-history  solution 


X  Vector  of  dependent  coordinates  resolved  to 

c 

global  reference  axis,  ND  x  1 


X  Coefficient  vector  of  first  sine  harmonic  term  in 

s 

harmonic  expansion  of  Xj  vector  (see  equation 
(24)),  N  x  1 


X  Coefficient  vector  of  sine  harmonic  term  in 

so 

harmonic  expansion  of  Xj  vector  for  no  HHC 
control  input  (see  equations  (84)  and  (85)), 
N  x  1 
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Y  Vector  of  coordinates  used  to  obtain  the  complex 

eigensolution  system  response  ( see  equation 
(35  )  ) ,  2N  x  1 

Vector  of  response  variables  in  the  expression 
for  the  vibration  reduction  performance  index 
(see  equations  (90)  and  (91)) 


Impedance  matrix  relating  rotor  hub  loads  to  hub 
displacements  (see  equation  (60)),  12  x  12 


Z11H'  Z12H' 
Z21H'  Z22H 


Submatrices  of  matrix  z„  (see  equation  (61)), 

ri 

6x6 


211D'  212D'  Submatrices  of  the  system  impedance  matrix 

2  ,  2  referred  to  dependent  coordinates  (see  equation 

r,  _L  L/  r,  /,  L/ 

(64)),  Nd  X  Nd 


P 


Transformation  matrix  relating  dependent  coordi¬ 
nates  to  independent  coordinates  (see  equation 
(4)),  ND  x  N 


3  Coefficient  of  cosine  term  in  harmonic  expansion 

c 

of  bifilar  mass  flap  displacement  (see  equation 
(77)),  rad 


6n  Flap  displacement  of  n-th  bifilar  mass  (see 

equation  (77)),  rad 

3  Coefficient  of  sine  term  in  harmonic  expansion  of 

s 

bifilar  mass  flap  displacement  (see  equation 
(77)),  rad 
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Transformation  matrix  relating  dependent  coordi¬ 
nates  resolved  to  substructure  local  orienta¬ 
tion  axis  to  dependent  coordinates  XQ  referred  to 
global  orientation  axis  (see  equation  (12)), 

ND  x  nd 

Bifilar  mass  flap  collective  displacement,  rad 

Transformation  ma  .rix  relating  coordinate  vector 
Xq  to  coordinate  vector  X^  (see  equation  (14)) 

Transformation  relating  coordinate  vector  to 

coordinate  vector  XJ2  (see  equation  (15)) 

Newmark-Beta  factor 


Modal  matrix  relating  displacements  at  a  node  to 
vector  of  normal  mode  coordinates  (see  equation 
(17)),  6  X  M 


Coefficient  of  cosine  term  in  harmonic  expansion 
of  bifilar  mass  lead  displacement  (see  equation 
(74)),  rad 


Lead  displacement  of  n-th  bifilar  mass  (see 
equation  (74)),  rad 


Coefficient  of  sine  term  in  harmonic  expansion  of 
bifilar  mass  lead  displacement  (see  equation 
(74)),  rad 


Fixed  absorber  mass  displacement,  m  (ft) 


e 


Rotor  blade  collective  pitch  angle  (scalar), 
rad 


e 

c 


e 


s 


<Vc 


<Vs 


Euler  angle  relating  local  and  global  axes  (see 
figure  ( 3 ) ) 

Vector  of  cosine  coefficients  of  HHC  control 
angles  (see  equation  (78)),  3  x  1 

Vector  of  sine  coefficients  of  HHC  control  angles 
(see  equation  (79),  3  x  1 

Coefficient  of  N-th  cosine  harmonic  in  blade  HHC 
control  angle  expansion 

Coefficient  of  N-th  sine  harmonic  in  blade  HHC 
control  angle  expansion 


91'  6 2'  03  Angular  displacements  of  a  connection  node  (see 

Figure  2) 


6,  <j> ,  1 1>  Euler  angles  relating  local  and  global  axes  (see 

Figure  ( 3  )  ) 

w  Frequency  of  a  normal  mode,  rad/sec 

Load  excitation  frequency,  rad/s 
0  Rotor  speed,  rad/s 
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Subscripts 


Superscripts 


Coefficient  of  cosine  term  in  a  harmonic  expan¬ 
sion  of  a  variable 


n-th  harmonic  of  stress  or  pushrod  load 


n-th  bifilar  mass 


time  index 


coefficient  of  sine  term  in  a  harmonic  expansion 
of  a  variable 


i-th  substructure 


imaginary  part 


index  of  last  substructure  in  a  set  of 
substructures 


real  part 


transpose 


inverse 


first  derivative  with  respect  to  time 
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second  derivative  with  respect  to  time 
amplitude 

optimum  value  of  a  control  angle 


OVERVIEW  OF  THE  COUPLED  ROTOR/AIRFRAME  VIBRATION  ANALYSIS 


The-  Coupled  Rotor/Airframe  Vibration  Analysis  consists  of  a  Base 
Program  which  is  used  in  conjunction  with  a  set  of  Externa] 
Programs  to  achieve  a  vibration  prediction  capability  for  heli¬ 
copters.  A  substructure  assembly  method  is  used  in  the  Base 
Program  to  assemble  a  dynamical  model  of  the  helicopter  from 
physical  helicopter  components  contained  in  the  Base  Program  and 
External  Program  models  of  the  rotor  system.  The  Base  Program 
solution  modes  include  steady  state  forced  responses,  eigen- 
solutions,  time  history  solutions,  and  parameter  variation  capa¬ 
bility.  Rotor  data  are  transferred  to  the  Base  Program,  which 
functions  in  a  stand-alone  mode  to  facilitate  the  evaluation  of 
parameter  variations  on  vibrations.  Stand-alone  operation  is 
achieved  with  a  small  core  requirement.  The  Base  Progiam  ana  che 
External  Programs  are  configured  to  provide  the  iser  with  latitude 
to  select  applications  tailored  to  user  inf- ,  aptitude,  and 
convenience.  If  rotor  data  are  on  file,  the  use  of  the  External 
Programs  is  not  required.  The  Base  Program  may  be  used  in  a 
stand-alone  mode  also  for  cases  not  requiring  rotor  data  -  e.g., 
investigation  of  transmission/isolator/fuselage  response.  Input 
to  the  Base  Program  is  designed  to  provide  an  economical  format 
requiring  input  only  for  selected  components.  Input  to  the  Base 
Program  is  not  restricted  to  the  External  Programs  discussed  in 
this  report.  Programs  with  similar  capabilities  may  be  utilized 
provided  that  input  formats  are  followed. 

Base  Program 

The  Base  Program  calculates  the  mass,  damping,  and  stiffness 
matrices,  and  forcing  function  vectors  of  substructures  resident 
in  the  Base  Program  code.  Input  data  consist  of  3-cha;acter 
literal  strings  which  define  the  components  to  be  included  in  the 
configuration.  These  3-character  literal  strings  trigger  input 
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requirements  for  the  components  and  the  forming  of  the  properties 
of  the  requested  components .  Input  data  also  define  component 
element  numbers  and  points  of  connection  with  other  components. 
The  Base  Program  uses  these  data  to  calculate  a  transformation 
matrix  which  synthesizes  the  coupled  equations  of  motion  for  the 
assembly  from  the  separate  equations  for  the  components.  The 
transformation  matrix  relates  the  set  of  coordinates  for  the 
components  of  the  requested  system  to  a  reduced  coordinate  set 
applicable  to  the  coupled  equations  of  motion.  Redundant  co¬ 
ordinates  are  removed  from  the  coordinates  for  the  requested 
components  (the  dependent  set)  to  create  the  coordinates  for  the 
coupled  system  (the  independent  set).  Rotation  transformations 
relating  component  local  axes  to  assembled  system  global  axes  are 
included  in  the  transformation  matiix.  The  user  may  elect  to 
solve  the  final  equation  of  motion  by  one  of  several  methods. 
Component  responses  are  recovered  from  the  solution  by  multi¬ 
plication  of  the  solution  vector  by  the  system  transformation 
matrix. 

The  External  Programs  yield  the  dynamic  and  aerodynamic  mass, 
damping,  and  stiffness  matrices,  and  forcing  functions  of  the 
rotor  system.  These  properties  are  assembled  with  the  properties 
of  the  components  resident  in  the  Base  Program  by  using  the 
transformation  method  described  above. 

Figure  1  illustrates  the  data  transfer  between  the  Base  Program 
and  the  External  Programs.  The  arrows  indicate  the  directions  of 
the  data  transfer.  Cycling  of  data  between  modules  is  indicated  by 
two  separate  arrows  joining  program  blocks.  The  character  of  the 
data  transferred  is  shown  in  Figure  1.  The  data  transfer  from  the 
External  Programs  into  the  Base  Program  is  accomplished  by  file 
transfers  and  the  Base  Program  operates  in  a  stand-alone  mode  to 
obtain  solutions. 
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Figure  1.  Coupled  Rotor/Airf rame  Vibration  Analysis  Data  Flows. 


The  Four  External  Programs  chosen  to  support  the  Base  Program  to 
form  the  coupled  rotor/airframe  analysis  are  described  in  a 
subsequent  section. 


SUBSTRUCTURE  ASSEMBLY  METHOD 

The  configuration  of  the  Base  Program  and  the  data  transfer 
between  the  Base  Program  and  the  External  Programs  are  based  on  a 
substructure  approach.  The  substructure  approach  synthesizes  the 
equations  of  motion  from  dynamical  models  of  parts  of  the  system. 
The  substructure  approach  is  the  key  to  a  flexible  modeling 
capability  permitting  wide  latitude  in  design  studies.  The 
approach  achieves  an  accurate  representation  of  rotor/airframe 
systems  in  the  framework  of  an  efficient  computer  program. 

The  substructure  method  employed  is  the  Hurty  method  of  Reference 
4.  The  coordinates  or  degrees  of  freedom  of  a  substructure  (also 
called  physical  component)  are  physical  and  generalized  displace¬ 
ments,  such  as  modal  amplitudes.  The  equation  of  motion  for  each 
substructure  is  expressed  in  mass,  damping,  and  stiffness  matrix, 
and  force  vector  form.  Character  and  treatment  of  exceptions  to 
these  forms  occurring  in  the  Base  Program  and  the  External  Pro¬ 
grams  are  described  subsequently.  Properties  of  the  i-th  sub¬ 
structure  are  [ml^^,  [c]^1^,  lk]^^,  {f}^^,  and  {rj  ^  ^ .  These 

are  respectively  mass,  damping,  and  stiffness  matrices,  and 
external  force  and  reaction  force  vectors.  The  submatrices  for 
the  substructures  are  collected  into  a  partitioned  diagonal  matrix 
equation  which  represents  the  system.  This  partitioned  diagonal 
matrix  equation  is 

Vd  +  CD*D  ♦  Vd  =  fd  -  "d  (1> 

The  matrix  M^  illustrates  the  typical  form  of  the  diagonal 
matrices  as: 
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(2) 


The  vector  of  coordinates  illustrates  the  form  of  the  vectors 
and  Rp  in  equation  (1). 


xD  =  !x!(l) 


The  subscript  D  signifies  that  equation  (1)  applies  to  a  set  of 
coordinates  ^  which  may  contain  coordinate  elements  which  are 
equal  to  or  depend  on  each  other.  In  other  words,  Xp  is  a 
"dependent"  coordinate  set,  which  can  contain  redundant  elements. 
The  coordinate  vector  Xp  contains  unconnected  coordinates  and 
connection  node  coordinates.  The  displacement  of  the  mass  of  a 
simple  absorber  would  be  an  unconnected  coordinate.  The  dis¬ 
placement  of  the  base  of  the  absorber  would  be  a  connection  node 
coordinate  because  the  base  of  the  absorber  may  be  attached  to 
another  component.  When  components  are  connected  to  each  other 
redundant  coordinates  occur  in  . 


A  force  vector,  {  f }  '  ,  represents  external  forces  and  moments 
applied  to  the  i-th  substructure.  The  forces  may  act  at  points  of 
the  structure  which  are  unconnected,  like  an  absorber  mass,  or  at 
a  connection  node,  like  the  absorber  base.  Reaction  forces  act  on 
a  substructure  at  points  of  connection  betwen  it  and  other  sub¬ 
structures,  and  are  induced  by  the  other  substructures.  The 
reaction  force  vector,  {r}^1^,  does  not  have  any  elements  asso¬ 
ciated  with  unconnected  coordinates.  Figure  2  shows  the  sense  of 
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(EXTERNAL  APPLIED  FORCE) 


Figure  2.  Displacements  and  Forces  Acting  on  a 
Connection  Node  of  a  substructure. 
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forces  and  moments  applied  to  a  connection  node.  Appendix  A 
illustrates  the  terminology  introduced  for  a  simple  system. 

The  synthesis  of  the  equations  of  motion  for  the  coupled  system  is 
accomplished  by  a  mapping  or  transformation  relating  the  dependent, 
coordinates,  X^,  to  a  reduced  coordinate  set,  .  Redundant 
coordinates  are  eliminated  by  requiring  component  displacements  to 
be  equal  at  connections .  The  reduced  coordinate  set  may  include 
modal  coordinates.  The  transformation  matrix  relating  X^  and  X^ 
is  denoted  by  p .  The  mapping  from  one  coordinate  set  to  the  other 
is : 


*D  =  P  XI 


Appendix  A  shows  a  sample  set  of  coordinates,  Xj^  and  Xj,  and  a 
sample  mapping. 

Substitution  of  the  transformation  equation  (4)  in  equation  (1) 

T 

and  pre-multiplication  of  equation  (1)  by  p  yields  the  final 
equation  of  motion  for  the  coupled  system.  This  is 


MXj  +  CX:  +  KXj.  =  F  (5) 

The  matrices  and  vector  in  equation  (5)  are 


PTMdP 

(6) 

pTcdp 

(7) 

ptV 

(8) 

^FD 

(9) 
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Matrix  p  is  assumed  to  be  a  constant  function  of  time,  since  its 
elements  relate  displacements  at  interfaces  which  are  rigidly 
joined,  or  coordinates  which  are  identical  -  e.g.  internal  co¬ 
ordinates.  The  condition  that  reaction  forces  at  connection  nodes 
are  equal  and  opposite  was  used  to  obtain  equation  (5).  This 
condition  reduces  to 


P^  =  0  (10) 

Equation  (10)  satisfies  the  requirement  that  work  done  by  internal 
reactions  of  the  assembled  structure  in  a  virtual  displacement, 
consistent  with  the  constraints,  is  zero. 

Equations  (6)  through  (9)  are  used  in  the  Base  Program  to  synthe¬ 
size  the  equations  of  motion  for-  the  coupled  system.  Appendix  A 
illustrates  the  synthesis  of  coupled  system  equations. 

The  component  displacements  X^  may  be  recovered  from  the  solution 
to  equation  (5),  Xj,  by  mapping  equation  (4).  Component  veloc¬ 
ities  and  accelerations  are  derived  similarly.  with  X^  and  its 
derivatives  known,  connection  node  or  interface  reactions  from 
equation  (1)  may  be  expressed  as 


*d  =  FD  -<mdxd  +  cd*d  +  W  (11> 

Rotor  hub  shears  and  moments  are  sample  interface  reactions. 
Appendix  A  illustrates  the  recovery  of  interface  reactions. 

Composition  of  Transformation  p 

The  transformation  matrix  p  is  expressed  as  the  product  of  three 
transformations  which  assemble  coupled  systems  from  components 
whose  properties  are  defined  in  component  local  axes,  and  allow 
for  the  use  of  modal  coordinates.  These  transformations  are 
described  in  the  following  subsections. 
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Transformation  from  Local  to  Global  Axes, 

The  transformation  from  dependent  coordinates  resolved  to  local 
axes  with  arbitrary  angular  orientations  to  dependent  coordinates 
resolved  to  a  global  reference  axis  is 


*D  =  PoXo  (12) 

Dependent  coordinate  sets  X^  and  Xq  can  include  coordinates  which 
are  redundant.  Redundant  coordinates  occur  when  components  are 
connected  to  each  other  at  connection  nodes.  Equation  (12) 
resolves  coordinates  at  connection  nodes  to  the  same  system  and  is 
a  prerequisite  to  the  elimination  of  redundant  coordinates  from 
X  .  Figure  3  illustrates  the  angular  orientation  of  a  loca!  axis 
system  relative  to  a  global  system.  The  local  system  is  obtained 
by  three  successive  rotation  displacements  of  axes  through  Euler 
angles  0  ,  $  ,  and  ij; . 

The  rotation  transformation  relating  linear  displacements  in  the 
local  and  global  axes  is 


Matrix  fE]  is  a  3  X  3  matrix  of  direction  cosines  each  of  which  is 
a  function  of  Euler  angles  0,  <|> ,  ij> .  The  same  transformations 
apply  also  when  small  angular  displacements  0^,  0^,  are  sub¬ 
stituted  for  u,  v,  w.  Linear  and  angular  velocities  and  ac¬ 
celerations  are  mapped  by  equivalents  of  equation  (13).  Figure  4 
defines  the  elements  of  matrix  [Ej.  Appendix  A  shows  the  ap¬ 
plication  of  a  rotation  transformation  to  a  simple  absorber.  The 
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absorber  mass  displacement,  6,  is  an  unconnected  coordinate  and  is 
retained  as  a  coordinate  in  the  final  system.  It  is  not  necessary 
to  resolve  6  to  a  global  system  because  unconnected  coordinates 
do  not  have  to  be  eliminated  by  equating  displacements  at  connec¬ 
tions.  Unconnected  coordinates  are  assumed  to  be  invariant  to 
rotation  transformations  from  local  to  global  axes. 

Transformation  from  Dependent  Coordinates  to  Physical  Domain 
Independent  Coordinates,  0^ 


The  transformation  from  dependent  coordinates  referred  to  a  global 
reference  axis  to  coordinates  from  which  redundant  coordinates  at 
connection  nodes  have  been  removed  is 

XQ  =  Px  XX1  (14) 

Coordinate  set  X^,  is  referred  to  as  the  physical  domain  inde¬ 
pendent  coordinate  set  since  connection  node  displacements  are 
physical  displacements.  The  determination  of  0^  is  based  on  a 
convention  in  the  computer  program  specifying  the  order  of  co¬ 
ordinates  in  the  dependent  sets  X^  and  X^.  Coordinates  in  Xq  are 
ordered  according  to  ascending  component  element  numbers  assigned 
by  the  user  in  the  input  specification  to  the  Base  Program. 
Appendix  A  illustrates  this  organization  in  XQ.  In  coordinate  set 
X j ^  ,  unconnected  coordinates  are  listed  first  in  order  of  ascend¬ 
ing  element  number.  Connected  coordinates  are  listed  after 
unconnected  coordinates  in  order  of  ascending  connection  mode 
numbers.  Redundant  coordinates  are  eliminated  by  first  retaining 
coordinates  for  a  component  with  the  lowest  element  number,  and 
then  eliminating  redundant  connected  coordinates  with  higher 
element  numbers  from  a  connection  node.  Appendix  A  shows  the 
ordering  of  coordinates  in  X^  .  The  logic  yielding  the  elements 

of  8 ,  is  based  on  the  ordering  of  coordinates  in  X  and  XT  , 
K 1  oil 
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described  above.  Appendix  A  shows  the  transformation  matrix  p^ 
corresponding  to  the  illustrations  given  for  XQ  and  . 

Transformation  from  Physical  Domain  Independent  Coordinates  to 
Modal  Domain  Independent  Coordinates, 


The  transformation  from  physical  domain  independent  coordinates, 

,  to  coordinates  which  include  modal  coordinates,  X^  /  is 


X1I  ^2  XI2 

This  transformation  is  required  when  modal  structures  are  used. 
Since  X j 2  ls  the  final  coordinate  set  for  the  system  to  be  solved 
it  is  symbolized  also  by  Xj ,  i.e. 


X 


I 


kI2 


(16) 


The  coordinate  set  XJ2  is  obtained  from  X^  by  expressing 
connection  node  coordinates  in  terms  of  modal  coordinates.  For 
example,  rotor/airframe  hub  coordinates  may  be  expressed  in  terms 
of  airframe  modal  coordinates.  Coordinate  set  X^2  may  contain 
connection  coordinates  which  have  not  been  expressed  as  modal 
coordinates,  and  will  in  general  include  unconnected  coordinates 
for  the  substructures.  Coordinate  set  XJ2  is  referred  to  as  the 
independent  coordinate  set  in  the  modal  domain. 

Connection  node  displacements  are  expressed  in  terms  of  modal 
coordinates  { q }  .  The  linear  and  angular  displacements  at  a  con¬ 
nection  node  referred  to  a  local  axis  for  the  modal  structure  are 
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u 


exl  =  ( Y  ]  {q! 


Local 


In  equation  (17), [y]  is  a  6XM  modal  matrix.  The  vector  of  normal 
modes  of  the  substructure  has  dimensions  MX1 ,  i.e.,  the  sub¬ 
structure  equations  of  motion  have  M  normal  modes.  Small  angular 
displacements  are  assumed.  Velocities  and  accelerations  are 
mapped  by  the  equivalents  of  equation  (17).  The  local  axis  system 
displacements  for  the  modal  structure  are  related  to  the  global 
axis  displacements  for  the  final  system  by 


[E]  ' 


- 1 - 


I  IEJ  I  O 


Local 


Global 


In  equation  (18),  [E]  is  the  matrix  of  direction  cosines  introduc¬ 
ed  in  equation  (13)  with  elements  defined  in  Figure  4.  To  obtain 
the  connection  node  displacements  referred  to  the  global  system  in 
terms  of  modal  coordinates,  the  inverse  of  the  rotation  transforma¬ 
tion  [E]  is  expressed  as: 


=  [E]  • 


?9 


Figure  4.  Elements  of  Rotation  Transformation  Matrix  [E]. 
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The  solution  to  equations  (17)  and  (18)  is 


T1 

[E]  .  0 

1 

l 


O 


Global 


i  T 

[El1 


[  Y  ]  { <3 } 


(20) 


Combining  equation  (20)  with  equation  (15)  yields  the  matrix  p^ 
relating  the  coordinates  in  to  modal  coordinates  and  other 
coordinates  in  X^.  Equation  (20)  is  used  at  each  connection  node 
where  physical  displacements  are  to  be  expressed  in  terms  of  modal 
displacements.  In  general,  the  matrices  [E]  and  [y]  will  differ 
among  the  nodes  because  different  modal  substructures  may  be 
joined  to  different  nodes  of  the  system,  and  because  of  different 
component  orientations.  Appendix  A  gives  examples  of  X^,  XJ2  and 

P2* 


Composite  Transformation  p  Relating  Dependent  Coordinates  to 
Independent  Coordinates 

From  equations  (12),  (14),  (15),  and  (16)  it  is  found  that 


^^o^l  P2  XI 


Comparison  of  equations  (4)  and  (21)  yields 


P 


P~  Pi  P- 


(22) 


31 


This  transformation  relates  the  dependent  coordinate  set,  X^,  for 
the  assembly  of  substructures  to  the  final  independent  (non- 
redundant)  coordinate  set,  Xj  ,  used  to  solve  the  equations  of 
motion.  Equation  (22)  is  used  to  construct  the  ^-transformation 
in  the  computer  program.  Appendix  A  gives  an  example  of  the 
composite  transformation. 


SOLUTION  MODES 


Four  solution  modes  of  equation  (5)  can  be  obtained.  Each  method 
has  different  attributes  affecting  information  provided,  speed  of 
execution,  and  computer  program  core  requirements.  These  solution 
modes  are:  1)  steady  state  forced  response,  FR1 ,  2)  eigensolution 
mode  for  systems  without  damping,  EG1,  3)  eigensolution  mode  for 
systems  with  damping,  EG2 ,  and  4)  time  history  mode  yielding 
system  response  as  a  function  of  time,  TH1 . 

In  addition  to  the  solution  modes  described  above,  two  more  opera¬ 
tional  modes  of  the  program  may  be  invoked.  These  modes  are  1) 
parameter  variation  mode,  PV1 ,  designed  to  provide  a  convenient 
selection  and  variation  of  parameters  in  design  studies,  and  2) 
general  program  control  mode,  GEN,  presently  configured  to  control 
display  of  program  output  and  graphics  options. 

Table  2  lists  the  character  names  associated  with  each  solution 
mode  and  operational  mode  available. 


A  description  of  each  solution  and  operational  mode  is  given 


TABLE  2  -  BASE  PROGRAM  SOLUTION  MODES  AND  OPERATIONAL  MODES 


DESCRIPTION  OF  MODE 


FORCED  RESPONSE  1  -  Forced  response  solution  mode 
yielding  harmonics  of  steady  state  response 

EIGENSOLUTION  1  -  Eigensolution  mode  for  systems 
without  damping 

EIGENSOLUTION  2  -  Eigensolution  mode  for  systems 
with  damping 

TIME  HISTORY  1  -  Time  history  solution  mode  using 
Newmark-Beta  method 

PARAMETER  VARIATION  1  -  Operational  mode  facilitat¬ 
ing  variation  of  parameters  for  design  studies 

GENERAL  -  General  control  mode  providing  general 
control  for  a  case 
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Steady  State  Forced  Response  ( FR1 ) 


The  harmonic  steady-state  forced  response  mode  FR1,  yields 
harmonic  coefficients  of  response  to  a  steady  harmonic  excit¬ 
ation.  This  solution  mode  is  useful  for  design  studies  by  virtue 
of  execution  time,  small  system  core  requirements,  and  importance 
of  steady  state  vibrations . 

The  forcing  function  which  appears  in  equation  (5)  is  assumed  to 
be 


F  =  F  cos  iuft  +  F  sin  ui-t  (23) 

CIS! 

In  equation  (23),  is  the  forcing  frequency  and  Fc  and  Fg  are 

N  X  1  vectors,  where  N  is  the  dimension  of  vectors  in  equation 
(5).  The  harmonic  steady  state  forced  response  is  also  assumed 
periodic . 


XT  =  X  cos  ui_t  +  X  sin  w,-t  (24) 
I  c  f  s  f  ' 

Substitution  of  equations  (23)  and  (24)  into  (5)  and  equating  the 
coefficient  expressions  of  cos  w^t  and  sin  w^t  to  zero  yields  the 
2N  X  2N  algebraic  equation 
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The  solution  to  equation  (25)  is 
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Determination  of  vibratory  response  requires  one  inversion, 
equation  (26).  Solution  speed  is  much  higher  than  a  time  history 
solution  and  the  method  is  well  suited  to  parametric  design 
studies . 

Equations  (26)  and  (26)  apply  to  systems  expressible  in  mass, 
damping,  and  stiffness  form.  These  matrices  may  contain  aero¬ 
dynamic  contributions  such  as  in  the  matrices  derived  from  aero- 
elastic  analysis  E927.  In  aeroelastic  analysis  G400,  rotor 
characteristics,  including  aerodynamic  effects,  are  represented  as 
a  hub  force  and  moment  excitation  vector  and  a  hub  impedance 
matrix.  Rotor  induced  impedance  matrices  must  be  added  to  the 
system  matrices  in  equation  (26)  to  obtain  the  vibratory  res 
ponse.  ihe  steady  state  solution  is 


-wf2M+K+Z1;L  J 


wfC+Z12 
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I  "wf  M+K+Z22 
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Forcing  function  vectors  Fc  and  Fg  in  equation  (27)  contain  rotor 
contributions.  Impedance  matrices  in  equation  (27)  apply  to  the 
assembled  system  with  independent  coordinates,  Xc  and  X  .  Trans¬ 
formations  similar  to  those  used  to  transform  mass,  damping,  and 

stiffness  matrices  in  equations  (6)  through  (8)  yield  the  imped¬ 

ance  matrices  in  equation  (27).  The  subject  of  mapping  rotor 
impedance  matrices  to  independent  coordinates  is  treated  in  a 
subsequent  section. 

The  system  vibratory  response  can  be  expressed  in  units  of  accel¬ 
eration  of  gravity  by 

xc/g  =  -(u»f2/9)Xc  (28) 

Xs/g  =  -(wf2/g)Xs  (29) 

The  computer  program  displays  the  amplitude  and  phase  of  each 
coordinate  in  X^  .  This  information  is  extracted  from  the 
elements  of  X^  and  X^  .  The  response  of  each  coordinate  is  assumed 
to  be  a  cosine  function.  The  standard  line-printer  output  does 
not  display  Xc  and  Xs>  A  debug  option  may  be  flagged  to  obtain  Xc 
and  Xs .  Amplitude  and  phase  data  only  are  stored  for  plotting. 

Eigensolution  for  Systems  Without  Damping  ( EG1 ) 

The  eigensolution  mode  for  systems  without  damping  yields  the 
frequencies  and  mode  shapes  of  structures  without  damping. 
Matrices  must  be  constant  with  respect  to  time.  The  effects  of 
changes  in  physical  properties  on  mode  shapes  and  frequencies  can 
be  determined.  These  effects  are  used  to  aid  the  interpretation 
of  the  effects  of  design  variables  on  normal  modes  and  vibratory 
response.  This  solution  mode  has  a  potential  capability,  not 
implemented  in  the  present  version  of  the  base  program,  to  trans¬ 
form  substructure  mass,  damping,  stiffness  and  force  vector 
properties  to  normal  mode  coordinates.  Normal  mode  coordinates 
reduce  the  problem  to  one  requiring  fewer  coordinates. 
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(30) 


Equation  (5)  is  specialized  to 
MX  +  KX  =  0 

The  solution  is  assumed  to  be  periodic 
Xj  =  Xelujt 

Substitution  of  equation  (31)  into  (30)  yields  the  NxN  algebraic 
equation 


(-w2M  +  K)  X  =  0  (32) 

Solutions  to  equation  (32)  yield  the  eigenvalues,  iu2,  and  cor¬ 
responding  eigenvectors,  X. 

Eigensoluticns  for  Systems  with  Damping  ( EG2 ) 


The  eigensolution  mode  for  systems  with  damping  is  used  to  provide 
information  on  the  stability  of  coupled  systems.  This  solution 
mode  is  restricted  to  hover.  Coupled  systems  may  include  rotor 
representations  which  embody  dynamic  and  aerodynamic  matrices 
derived  from  rotor  aeroelastic  analysis  E927.  Stability  informa¬ 
tion  may  be  used  to  supplement  the  steady  state  vibratory  response 
data  for  a  system  to  provide  more  data  on  system  attributes. 
Matrices  must  be  constant  with  respect  to  time. 

Equation  (5)  is  specialized  to 

MXj  +  CXj  +  KXj  =  0  (33) 

Reduced  to  the  first  order  system  equation  (33)  becomes 
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AY  -  BY  =  O 

This  is  accomplished  by  the  substitution 


Matrices  in  equation  (34)  are 
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The  eigenvalue  problem  is  obtained  from  equation  (34)  by  assuming 


Y  =  Ye  At  (38) 

which  yields  the  2N  X  2N  algebraic  equation 

(\  A  -  B ) Y  =  O  (39) 

The  modal  frequency  in  hertz  and  modal  damping  in  fractions  of 
critical  damping  respectively  are 


with 


w  =  \^/(2n) 

C  =  +  Ar2,H 


\  =  xr  +  i  Ki 


(40) 

(41) 


(42) 
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representing  a  complex  eigenvalue.  The  real  and  imaginary  com¬ 
ponents  of  Y  are  displayed  for  each  mode.  A  stable  mode  corres¬ 
ponds  to  a  positive  value  of  £  (negative  value  of  A  ). 

Time  History  Solution  Mode  (TH1) 

The  time  history  solution  mode  TH1 ,  yields  the  response  of  a 
system  as  a  function  of  time.  Transient  responses,  steady 
responses,  and  the  stability  of  systems  can  be  determined  from  the 
time-history  behavior  of  a  system.  Responses  of  nonlinear  systems 
can  be  determined.  The  time-history  mode  is  able  to,  or  has  the 
potential  to,  supply  all  the  information  provided  by  forced 
response  and  eigensolution  modes,  in  addition  to  the  analysis  of 
transient  and  non-linear  data.  A  drawback  of  the  method  is  its 
high  execution  time  in  comparison  with  the  other  methods  making  it 
comparatively  costly  and  inefficient  for  design  studies. 

The  solution  algorithm  yielding  the  time  history  response  is  the 
Newmark-Beta  method  described  in  Reference  5.  The  method  is  a 
finite  difference  method.  Displacement  responses  are  obtained 
from  displacements  known  at  prior  times,  and  from  data  defining 
the  magnitudes  of  known  external  forces  acting  on  the  system. 

The  vector  of  coordinates  satisfying  equation  (5)  at  time  step  n 
is  denoted  (X^ )n»  The  corresponding  time  is 


t 

n 


t  ,  +  nh 
n-1 


(43) 


In  equation  (43)  h  is  the  step  size.  The  integer  step  number  n 
ranges  from  one  to  the  number  of  steps  in  the  calculation.  Time 


t  is  the  initial  time. 

<Vo  and  <Vo- 


Corresponding  initial  conditions  are 
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When  (Xj)o  and  (xj)Q  are  specified  the  displacement  response  at  t 
is  found  from 

(XI)1  -  D_1  [P(XI)Q  +  Q(Xi)q  +  fS*h2Fx  +  RFq]  (44) 

Matrices  D,  P,  Q  and  R  are  functions  of  M,  C,  and  K  in  equation 

(5),  time  step  h,  and  factor  0*  (see  equations  (48)  through  (52)). 

They  have  dimensions  NXN  where  N  is  the  size  of  Xj  in  equation 

(5).  The  Factor  0*  is  the  Newmark-Beta  Factor.  Values  of  0* 

ranging  from  O  to  0.25  can  be  input  to  the  Base  Program.  The 

factor  0*  is  used  to  control  the  variation  of  acceleration  assumed 

in  the  time  interval  (see  remarks  in  Reference  5).  A  numerical 

solution  with  0*  =  0.25  is  unconditionally  stable.  Force  vectors 

F,  and  F  are  known  external  forcing  functions  at  t  and  t, . 

1  o  o  1 

After  the  solution  (X^J^  is  known,  successive  solutions  are 
obtained  from  the  recursion  formula 

<V„+1  =  0-1  '  (xl'n  -F'<xI>n-l 

+  0*h2  <Fn+1  +  (l/0*-2)Fn  +  Fn_1)]  (45) 

The  NxN  Matrices  B  and  F'  are  known  functions  of  M,  C,  and  K,  step 
size  h,  and  factor  p*  (see  equations  (53)  and  (54)). 

Equation  (45)  may  be  used  also  to  restart  solutions  from  solutions 
calculated  up  to  the  time  of  restart. 


The  additional  assumption  is  made  in  the  Base  Program  code  that 

the  force  F  , ,  linearly  extrapolates  the  forces  F  and  F  For 

n+1  n  n-1 

equal  intervals  of  time  h  this  assumption  reduces  to 


F 


n+1 


2Fn 


n-1 


(46) 
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and  equation  (44)  becomes 


-1 


<XI>n+l  =  D 


tB<XI>n  -  F'<XI>n-l  +  h2Fn 


(47) 


Equation  (46)  allows  for  the  inclusion  of  aerodynamic  forces  which 
are  functions  of  Xj  ,  and  .  The  present  version  of  the  Base 
Program  does  not  model  these  aerodynamic  forces.  A  technique  to 
include  these  forces  is  to  assume  F  to  be  a  function  of  Xj  and  its 
derivatives  at  a  prior  time.  Equations  (46)  and  (47)  are  consis¬ 


tent  with  this  assumption. 

Matrices  in  equations  (44),  (45),  and  (47)  are 

D  =  M  +  ( h/2 ) C  +  0*h2K  (48) 

P  =  M  +  ( h/2 ) C  -  ( 1/2-0 * )  h2K 

- ( 1/4-0  * )h3CM_1K  (49) 

Q  =  [M  -  (1/4  -  0*)  h2CM-1C ]  h  (50) 

R  =  [(1/2  -  0*)I  +  (1/4  -p*)  hCM-1  ]  h2  (51) 

I  =  unit  matrix  (52) 

B  =  2M  -  (1-20*)  h2K  (53) 

F'  =  M  -  ( h/2 ) C  +  0 *h2K  (54) 


The  dimensions  of  matrices  in  equations  (48)  through  (54)  are  NXN 
where  N  is  the  dimension  of  the  vector  of  independent  coordinates 

xr 
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Parameter  Variation  Mode  (PV1) 


The  parameter  variation  mode,  PV1 ,  is  used  to  identify  and  specify 
a  range  for  a  parameter  which  is  to  be  changed  in  design  studies. 
This  operational  mode  provides  a  convenient  and  general  mode  of 
selection  of  parameters  to  be  varied  among  the  substructures 
available  in  the  Base  Program.  Up  to  ten  parameters  from  any  of 
the  substructures  may  be  varied  simultaneously. 

General  Control  Mode,  GEN 

The  general  control  mode,  GEN,  is  designed  to  provide  a  basis  for 
the  application  of  general  controls  to  a  case  such  as  line  printer 
output  and  graphics  selection. 


ASSEMBLY  OF  ROTOR  IMPEDANCE  MATRIX 

The  purpose  of  this  section  is  to  derive  the  procedure  for  the 
assembly  of  a  rotor  impedance  matrix  with  mass,  damping,  and 
stiffness  matrices  and  forcing  functions  of  other  substructures. 
The  mapping  of  the  rotor  impedance  matrix  to  independent  coordi¬ 
nates  for  the  coupled  system  is  shown  to  be  the  same  as  for  other' 
substructure  properties .  The  rotor  impedance  matrix  must  be 
partitioned  to  accomplish  the  mapping. 

The  effects  of  hub  displacement  are  included  in  he  coupled 
rotor/airframe  system  by  representing  the  rotor  effects  as  the  sum 
of  a  hub  excitation  for  a  fixed  hub  and  a  correction  involving  the 
rotor  impedance  matrix  multiplied  by  hub  displacement.  This 
representation  enables  rotor  properties  to  be  calculated  prior  to 
execution  of  the  Base  Program.  This  formulation  maintains  the 
desired  small  core  and  execution  time  for  the  Base  Program. 
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To  obtain  a  rotor  hub  excitation  vector  and  a  rotor  hub  impedance 

matrix,  the  rotor  is  first  analyzed  for  a  fixed  hub  in  the  aero- 

elastic  analysis  -  e.g.  G400.  The  six  hub  forces  and  moments 
acting  on  the  fixed  hub  are 

US  H  =  UCS  H  cos  U)ft  +  {  f  s  5  H  sin  wrt  (55) 

Vectors  { f  }  and  {f  }  „  are  6X1  component  vectors  of  cosine 

C  tt  S  Jti 

and  sine  terms  in  equation  (55).  For  a  N-bladed  rotor  with 
identical  blades,  frequency  is  equal  to  nNO ,  where  n  is  an 

integer  and  0  is  the  rotor  speed.  Forces  and  moments  at  frequency 
are  the  only  vibratory  forces  that  can  be  transmitted  to  the 
non-rotating  axis  at  the  hub. 

To  calculate  elements  of  the  rotor  impedance  matrix  in  aeroelastic 
analysis  G400,  the  rotor  hub  is  perturbed  harmonically  in  the 
fore-aft  direction  by  a  displacement  uccos  w^t.  This  displacement 
causes  a  force  perturbation  acting  on  the  hub.  The  reaction  force 
acting  on  the  hub  can  be  expressed  as 

{r}  H  =  {fc}  H  COS  “f1  +  {fs]  H  Sin  Wffc 

■*zcl*  UC  COS  wft  “  *zsl*  UC  sin  wft  (56) 

The  6X1  reaction  force  is  equal  to  the  constrained  hub  forces 
augmented  by  terms  which  are  proportional  to  the  hub  fore-aft 
displacement,  uc*  The  6X1  column  vectors  ?zcll  and  izslJ  are 
columns  of  the  impedance  matrix.  Elements  in  these  vectors  are 
assumed  to  be  constant.  In  a  similar  way,  lateral  perturbations 
of  the  hub,  vccos  w^t,  are  introduced  to  derive  impedance  columns 
{z  }  and  {z  _}  .  This  procedure  is  repeated  for  all  possible 
displacements  at  the  hub,  including  displacements  proportional  to 
sin  w^-t.  The  complete  rotor  impedance  matrix  and  the  expression 
for  the  hub  reactions  including  all  hub  displacements  is  obtained. 
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The  following  definitions  express  harmonic  coefficients  for  the 
hub  reactions  as  functions  of  the  impedance  matrix.  The  6X1  hub 
reactions  are  assumed  to  be 

{r}  H  =  {rc}  H  cos  uift  +  { rs }  R  sin  wft  (57) 

The  6X1  hub  displacement  perturbation  vector  is 

{x}  H  =  {xc?  H  cos  wft  +  {xs}  H  sin  wft  (58) 


with 


T 


T 


{*}  H  =  {u,  v,  w,  e1(  &2,  e3) 


(59) 


Expressions  similar  to  equation  (59)  define  {x  )  „  and  {x  }  „ 

C  ri  S  H 


Equating  the  coefficient  expressions  for  cos  >u£.t  and  sin  u>  t 
in  the  reaction  force  equation  tc  zero  yields 
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j  Z22H 
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(60) 

(61) 


Matrix  z  is  a  12  X  12  rotor  hub  impedance  matrix  composed  of  four 

n 

6X6  submatrices  Z^2H'  Z21H'  Z22H  as  shown  equation 

(61).  The  hub  reaction  force  in  equation  (60)  is  composed  of 
fixed  hub  forces  {f  }  and  {f  }  „  and  a  correction  proportional 
to  the  impedance  matrix  z„  multiplied  by  the  hub  displacement 

ri 

vector.  Substituting  equation  (61)  into  (60)  and  rearranging  the 
terms  yields 
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To  develop  the  algorithm  for  assembling  impedance  properties  in 
equation  (62)  with  properties  of  other  substructures  expressed  in 
mass,  damping,  and  stiffness  form,  a  rotor  impedance  substructure 
is  introduced.  This  has  null  mass,  damping,  and  stiffness 
matrices,  and  non-zero  hub  force  and  moment  excitation  and  rotor 
hub  impedance.  The  vector  of  the  dependent  coordinates  X^  is 
expressed  as 

*D  =  <xD)c  COS  ^f*1  +  ^ XD ^ s  Sin  Lft  (63) 

External  forces  and  reactions  RD  are  expressed  similarly. 
Substructure  mass,  damping,  and  stiffness  matrices  are  collected 
with  impedances  to  yield 


--  w  f -  kd  +  z11D  | 

1 

WfCD  +  Z12D 

UJfCD  +  Z21D  j’ 

Wf2MD  +  ^  +  Z22D 
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.  \^J 

<FD>s( 

K>s\ 

Vectors  ( X  )  and  (X_)  contain  rotor  hub  coordinates  jx  }  tt  and 
Dc  us  c  H 

{x  i  „  from  equation  (62).  In  equation  (b4),  the  partitioned 

S  n 

matrices  are  diagonal.  Matrices  Mp,  Cp,  and  Kp  are  typified  by 
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md  = 
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(65) 


Zero  diagonal  submatrices  [Oh  are  associated  with  coordinates 

H 

{x  }  and  [x  }  which  have  no  mass,  damping,  or  stiffness 

C  H  S  H 

matrices  in  the  rotor  impedance  substructure  defined  by  equation 
(62).  The  impedance  matrix  Z  in  equation  (64)  is 


(66) 


The  matrix  in  equation  (66)  is  zero  except  for  a  6  X  6  diagonal 

element  z11u  associated  with  hub  vector  {x  Element  z 

11H  CH  11  rl 

derives  from  the  impedance  substructure  defined  in  equation  (62). 
The  forms  of  impedance  matrices  Z  12d'  z21D'  Z22D  111  e<3uat:on  (64) 
are  similar  to  equation  (66). 

Substructure  RE3  which  models  a  rotor  hub  excitation  force  and  a 
hub  impedance  matrix  incorporates  the  properties  described  above. 
Each  6X6  mass,  damping,  and  stiffness  matrix  for  RE3  is  zero. 
Substructure  impedance  elements  and  hub  excitation  forces  are 
transmitted  from  external  aeroelastic  analysis  G400. 

The  derivation  of  the  transformation  to  independent  coordinates  is 
accomplished  by  substituting  equation  (63)  for  X^  and  equation 
(24)  for  Xj  into  equation  (4).  This  yields 
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Substitution  of  equation  (67)  and  (68)  in  equation  (64)  and  pre¬ 
multiplication  of  equation  (64)  by 
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yields  equation  (27)  for  the  vibratory  response  of  the  coupled 
system  given  by 


i7\ 


are 


--  Wf* 

M+K+Z1:l 

1 

+  u*fC+Zl2  " 

F  . 

\  c  i 

(27) 

1 

-  ui  „C+Z„  . 

I 

-  w  <r2M+K+Z_ 

If  ' 

L  1 

21 

\ 

f  22 

s 

x  mappings 

yieldi 

ng 

the  matrices  in  equation 

(27) 

211 

= 

"XlD 

P 

(70) 

212 

= 

f>Tzi2D 

P 

(71) 

Z12 

= 

P 

(72) 

„T_ 

(  73  ) 

Z22 

— 

P  Z22D 

P 

47 


In  the  computer  program  implementation  the  zero  mass,  damping  and 
stiffness  properties  of  substructure  RE3  are  assembled  with  the- 
properties  of  other  substructures  to  form  Mp,  C^,  and  in 

equation  ( 1 ) .  These  matrices  are  mapped  to  the  independent 
coordinate  system  by  equations  (6)  through  (8).  The  forced 
response  solution  mode  is  invoked  and  matrices  in  equation  (26) 
are  formed.  Impedance  matrices  are  mapped  to  independent  co¬ 
ordinates  by  equations  (70)  through  (73)  and  are  added  to  the 

matrices  of  equation  (26)  to  form  equation  (27).  Equation  (27) 
yields  the  vibratory  response  of  systems  including  a  rotor  imped¬ 
ance  model . 

Hub  reactions  can  be  recovered  from  equations  (67),  (68),  and  (64) 

after  the  coupled  system  solution  is  found  from  equation  (27). 

The  RE3  impedance  substructure  is  consistent  only  with  a  steady 
state  forced  response  mode  of  solution,  FR1 .  The  forcing  fre¬ 
quency  in  FR1  must  be  the  same  as  the  forcing  frequency  \ised  to 
derive  rotor  impedance  properties.  This  frequency  will  be  an 

integer  multiple  of  rotor  speed  and  blade  number  (w^  =  nN  Q). 

BASE  PROGRAM  SUBSTRUCTURES 

Substructures  available  in  the  Base  Program  listed  in  Table  3  are 
described  below.  Descriptions  are  provided  of  substructure 
function,  coordinates,  properties,  and  assumptions. 

Substructure  FA1 

The  function  of  substructure  FA1  is  to  simulate  a  vibration 
absorber  which  may  be  attached  to  any  other  substructure  in  the 
fixed  -  i.e.,  non-rotating-system .  Substructure  FA1  has  seven 

coordinates.  These  are  absorber  mass  displacement,  and  six 
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TABLE  3  -  BASE  PROGRAM  SUBSTRUCTURES 

COMPONENT 

FLAG 

DESCRIPTION  OF  COMPONENT 

FA1 

FIXED  ABSORBER  1  -  Single  mass,  spring  and  damper 

vibration  absorber. 

GF1 

GENERALIZED  FORCE  1  -  Force  component  allowing  speci¬ 
fication  of  harmonic  force  and  moment  excitation. 

BM1 

BEAM  MODEL  1  -  Uniform  elastic  beam  segment. 

CN1 

CONSTRAINT  COMPONENT  1  -  Component  allowing  elimina¬ 
tion  of  coordinates  (degrees  of  freedom)  from  a 

connection  node. 

RE  2 

ROTOR  ELASTIC  2  -  Model  of  aeroelastic  rotor  expressed 
in  mass,  damping,  and  stiffness  form.  Rotor  data  are 

transmitted  to  the  Base  Program  by  rotor  aeroelastic 

analysis  E927.  Restricted  to  hover. 

MSI 

MODAL  SUBSTRUCTURE  1  -  Normal  mode  representation  of  a 

dynamical  structure. 

BF1 

BIFILAR  1  -  Horizontal  linear  bifilar  vibration 

absorber  for  reducing  rotor  inplane  vibrations. 

BF2 

BIFILAR  2  -  Vertical  linear  bifilar  vibration  absorber 

for  reducing  rotor  out  of  plane  vibrations. 

i  si 

ISOLATOR  1  -  Nodal  isolator  employing  two  anti -resonant 

bar  weights. 

RE3 

ROTOR  ELASTIC  3  -  Impedance  model  of  aeroelastic  rotor 

in  hover  or  forward  flight.  A  hub  force  and  moment 

excitation  vector  and  a  hub  impedance  matrix  are  trans¬ 
mitted  to  the  Base  Program  by  rotor  aeroelastic  analysis 

G400.  Vibratory  response  of  a  coupled  system  is  found 

using  solution  mode  FR1 .  This  component  has  a  cap¬ 
ability  to  obtain  optimal  HHC  control  angles  for  vibra¬ 
tion  reduction. 
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displacements  of  the  base  consisting  of  three  linear  base  dis¬ 
placements  and  three  angular  base  rotations.  The  mass  of  the 
absorber  is  free  to  move  perpendicular  to  the  base  and  does  not 
move  parallel  to  the  base.  Viscous  damping  is  assumed.  Figure  5 
illustrates  substructure  FA1  properties.  Initial  mass  displace¬ 
ment  and  velocity  can  be  specified  for  transient  response  ana¬ 
lysis  . 


Substructure  GF1 

Generalized  Force  1,  GF1 ,  specifies  a  harmonic  force  and  moment 
excitation  acting  at  a  point  on  a  structure.  The  GF1  substructure 
has  six  coordinates,  comprising  three  linear  displacements  and 
three  angular  rotations  at  a  connection  node.  Forces  and  moments 
are  resolved  to  a  non-rotating  axis  system.  The  user  specifies 
the  cosine  and  sine  coefficients  of  these  forces  and  moments  to 
define  a  harmonic  excitation  of  known  magnitude.  The  mass, 
damping,  and  stiffness  matrices  for  substructure  GF1  are  null 
matrices.  Substructure  GF1  may  be  used  to  represent  an  excitation 
for  forced  response  (FR1)  or  time-history  ( TH1 )  solution  modes. 
Substructure  GF1  may  be  included  in  eigensolution  (EG1  and  EG2 ) 
input  streams  without  affecting  the  final  eigensolution.  Figure  6 
illustrates  substructure  GF1  properties. 

Substructure  BM1 

Beam  Model  1,  BM1 ,  is  a  uniform  elastic  beam.  This  substructure 
has  12  coordinates  consisting  of  three  linear  displacement  and 
three  angular  rotations  at  each  end  of  the  beam.  The  beam  may  be 
connected  to  other  substructures.  The  uniform  mass  of  the  beam, 
mo,  is  approximated  by  two  concentrated  masses  equal  to  mQ/2  at 
its  ends.  The  beam  is  assumed  to  have  zero  mass  moments  of 
inertia  associated  with  bending  and  torsion  rotations.  Structural 
damping  is  assumed  to  be  zero.  Figure  7  illustrates  substructure 
BM1  properties. 
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Figure  7.  Uniform  Beam  Segment  BML 


Constraint  Substructure  CN1 


Constraint  substructure  CN1  may  be  used  to  eliminate  or  constrain 

coordinates  at  a  connection  node.  There  are  six  coordinates  at 

each  connection  node  (see  Figure  2).  The  coordinates  consist  of 

three  linear  displacements  and  three  angular  rotations.  The  user 

may  constrain  any  five  of  the  six  coordinates  at  a  connection 

node.  At  least  one  coordinate  among  the  six  must  be  retained. 

The  constraint  component  is  a  specialized  modal  substructure  (see 

Substructure  MSI )  whose  modal  matrix  contains  ones  and  zeroes . 

The  occurrence  of  a  one  in  a  modal  matrix  row  causes  a  connection 

node  coordinate  to  be  replaced  by  a  modal  coordinate  which  is 

identical  to  the  connection  coordinate.  The  occurrence  of  zeroes; 

in  a  modal  matrix  row  eliminates  a  connection  node  coordinate. 

Equation  (20)  explains  how  the  connection  node  coordinates!  are 

T 

eliminated.  Matrix  [E]  ,  containing  the  rotation  transformation 
from  local  to  global  axes  is  replaced  by  the  identity  mat  i  ix  I  . 
The  number  of  coordinates  in  vector  {a!  is  selected  to  he  equal  to 
the  retained  coordinates  at  the  connection  node.  The  ami  d ;  t.ut  es. 
in  vector  {q}  are  identical  to  coordinates  to  be  tetained  at  the 
connection . 

Rotor  Aeroelastic  Substructure  RE2  and  External  liogta.'i  !'•. 

Rotor  aeroelastic  substructure  RE2  is  used  to  represent  a  multi¬ 
blade  aeroelastic  rotor  in  hover.  This  component  identities  rotor 
matrices  transmitted  from  the  rotor  aeroelastic  analysis  K92'7  and 
arranges  the  rotor  data  for  assembly  with  ottiei  components. 
Standard  application  of  RE 2  is  in  conjunction  with  eigensolut ion 
mode  EG2 . 

Table  4  lists  RE2  coordinates,  expressed  as  output  labels  dis¬ 
played  by  the  Base  Program.  Use  of  rotor  coordinates  in  the  rotor 
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TABLE  4 


AEROELASTIC  ROTOR  RE2  COORDINATE  LABELS 


COORDINATE 

DESCRIPTION 

OUTPUT  LABEL 

QFT1-QFT4 

Blade  bending  collective  modal 

amplitude  ^  a  ^ 

THTR 

Blade  pitching^*5) 

THTE 

Blade  torsion^  ^ 

BETA 

/  c ) 

Flapping  collective  mode 

GAMA 

Leading  collective  mode^^ 

QT1 S ,  QT1C 

Blade  bending  cyclic  sine  and  cosine 

QT2S ,  QT2C 

modal  amplitudes'3^ 

QT3S,  QT3C 

QT4S ,  QT4C 

THRS ,  THRC 

Blade  pitching  sine  and  cosine  cyclic 

amplitudes^  ^ 

THES ,  THEC 

Blade  torsion  sine  and  cosine  cyclic 

amplitudes  ^  ^ 

BETS ,  BETC 

Blade  flapping  sine  and  cosine  cyclic 

( c ) 

amplitudes 

GAMS ,  GAMC 

Blade  leading  sine  and  cosine 

amplitudes  ^  ^  ^ 

X,  Y,  Z,  THTX , 

Hub  displacements.  See  Figure  2  for 

THTY,  THTZ 

positive  directions 

NOTES : a )  Positive  flatwise  bending  up  and  edgewise  bending 

forward . 

b)  Positive  nose  down. 

c)  Positive  up. 

d)  Positive  forward. 

equations,  and  the  assumption  of  identical  blades,  equally  dis¬ 
posed  on  the  rotor  disc,  enables  a  constant,  coefficient  equation 
set  to  be  derived  (see  Reference  6).  Rotor  coordinate  selection 
is  made  in  the  input  specification  to  E927.  Rotor  aeroelastic 
analysis  E927  transmits  to  a  file  the  number  of  coordinates,  flags 
defining  coordinates  selected,  and  corresponding  mass,  damping, 
and  stiffness  matrices.  Matrices  embody  dynamic  and  aerodynamic 
contributions.  Numerical  input  data  following  the  RE2  flag 
consist  of  hub  connection  node  number  and  hub  orientation  angles. 
Processing  by  the  Base  Program  includes  the  addition  of  a  dummy 
yaw  coordinate.  Zero  mass,  damping,  and  stiffness  matrices  are 
created  for  this  coordinate.  This  yields  a  six  coordinate  hub 
consistent  with  coordinate  representation  for  other  substructures. 
Six  hub  coordinates  enable  the  rotor  component  to  be  coupled  to 
other  components  without  changing  Base  Program  assembly  logic.  A 
valid  physical  model  of  rotor  yaw  motion  is  not  achieved.  Null 
rotor  matrices  are  associated  with  this  displacement.  The  user 
must  eliminate  the  hub  yaw  coordinate  from  the  system  model.  This 
may  be  accomplished  with  constraint  component  CN1 .  A  valid  model 
would  occur  also  if  fuselage  normal  modes  in  component  MSI  have  no 
hub  yaw  displacement. 


Modal  Substructure  MSI 

Modal  Substructure  MSI  represents  a  dynamical  system  employing 
normal  modes  of  a  substructure.  Each  substructure  of  MSI  type  is 
limited  to  the  use  of  one  normal  mode.  The  user  must  supply  input 
defining  the  modal  generalized  mass,  damping  ratio,  and  frequency 
of  the  normal  mode.  These  properties  can  be  generated  from 
analysis  or  from  test  data.  Processing  in  the  program  forms 
damping  and  stiffness  of  the  mode  from  the  input  data.  Each 
substructure  of  type  MSI  is  allowed  to  have  up  to  five  nodes  at 
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which  node  number,  structure  orientation  angles,  and  elements  of 
the  modal  matrix  for  the  normal  mode  are  defined.  At  each  such 
node,  the  elements  of  the  modal  matrix  are  three  linear  d^  ’  ace- 
ments  and  three  angular  rotations  for  a  unit  value  of  tht  .ormal 
mode  amplitude.  The  directions  of  the  displacements  at  a  node 
must  be  consistent  with  the  convention  shown  in  Figure  2.  The 
specification  of  the  same  node  number  between  two  substructures  of 
any  type  leads  to  the  connecting  of  the  two  substructures  at  the 
node.  More  than  one  modal  substructure  of  type  MSI  is  connected 
to  the  same  node  when  it  is  desired  to  represent  the  effects  on  a 
node  of  several  normal  modes.  Nodes  on  MSI  substructures  n:r  y  be 
used  to  monitor  airframe  responses  at  up  to  five  points.  If  a 
physical  substructure  is  not  connected  to  a  node  at  which  re¬ 
sponses  are  desired,  these  responses  may  be  displayed  by  con¬ 
necting  a  GF1  substructure  to  the  node  of  interest,  with  zero 
forces  and  moments  specified.  This  GF1  substructure  contributes  no 
matrices  or  forces  to  the  system.  It  introduces  physical  dis¬ 
placements  at  the  node  which  become  members  of  the  set  of  depend¬ 
ent  system  coordinates  which  are  used  to  monitor  the  node  response 
Initial  modal  amplitudes  and  velocities  can  be  specified  to 
analyze  system  transient  response. 

Horizontal  Linear  Bifilar  BF1 

Horizontal  linear  bifilar  BF1  represents  a  bifilar  vibration 
absorber  designed  to  reduce  N  per  rev  vibrations  in  the  plane  of 
rotation  of  a  helicopter  rotor.  The  bifilar  absorber  is  assumed 
to  consist  of  N  arms  equal  in  number  to  the  rotor  blades.  Each 
arm  holds  a  pendular  mass  free  to  lead  or  lag.  The  arms  are 
assumed  to  be  equally  disposed  and  each  holds  an  identical  ab¬ 
sorber  mass.  The  bifilar  absorber  rotates  at  the  helicopter  rotor 
speed  0.  Figure  8  is  a  schematic  of  one  arm  of  the  bifilar.  Each 
bifilar  mass  lead-lag  displacement  may  be  expressed  as  the  sum  of 
bifilar  lead-lag  collective  and  cyclic  displacements.  The  lead- 
lag  displacement  vn  of  the  n-th  bifilar  is 
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N  MASSES,  ONLY  ONE  MASS  SHOWN 
y  =Ky  COS<^  +  y  SIN  \ft  )(2/N) 

C  5 _ _ 


Figure  8.  Horizontal  Linear  Bifilar  BF1. 
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(74) 


Y n  -  yo/N  +  ( V ccostii^  +  Yssin  i\>n)(2/N) 

4>n  =  (2rx  /N)  ( n-1  )  +  «|>  (75) 

«(:  =  fit  (  76  ) 

Equation  (74)  introduces  the  bifilar  collective  and  cyclic  coordin¬ 
ates  in  the  equations  of  motion  of  the  bifilar.  The  angle  o  is 
the  azimuth  angle  of  the  n-th  bifilar  arm.  The  equations  of 
motion  for  each  bifilar  consist  of  N  equations  defining  the 
equilibrium  of  moments  about  the  hinge  point  of  each  pendular 
mass.  These  N  equations  are  supplemented  by  six  equations  defin¬ 
ing  the  equilibrium  of  the  bifilar  hub  in  three  translational  and 
three  rotational  directions.  When  equation  (74)  is  substituted 
into  these  equations,  an  equation  set  with  eight  coordinates  is 
obtained.  The  assumption  is  made  that  yQ  =  0.  The  eight  co¬ 
ordinates  consist  of  two  bifilar  mass  sine  (y  )  and  cosine  (y  ) 

s  c 

cyclic  amplitudes  and  six  hub  displacements  shown  in  Figure  8. 
The  reduced  equation  set  has  constant  coefficients.  Restrictions 
on  the  configurations  that  can  be  analyzed  are  bifilars  which  have 
three  or  more  equally  disposed  and  equal  masses.  Pendulum  dis¬ 
placements  are  assumed  to  be  small  enough  to  permit  linearizations 
of  the  equations.  Solutions  are  obtained  in  terms  of  the  cyclic 
amplitudes  vc  and  y  ,  as  functions  of  time.  Solution  modes  that 
may  'be  used  with  BF1  are  FR1 ,  EG1 ,  EG2 ,  and  TH1 .  Initial  values 
of  bifilar  mass  cyclic  mode  amplitudes  can  be  specified  for  tran¬ 
sient  analysis. 

Horizontal  bifilars  are  tuned  to  eliminate  hub  vibrations  induced 
by  forces  referred  to  the  rotating  system  at  frequencies  (N-1)  Q 
or  (N+l)  Q,  respectively.  These  forces  are  transformed  to  forces 
with  frequency  NQ  resolved  to  non-rotating  axes.  Thus,  two  separ¬ 
ate  horizontal  bifilars  are  required  to  eliminate  NQ  forces  re¬ 
solved  to  the  fixed  system.  Post-processing  of  response  results 
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obtained  from  the  system  solution  may  be  required  to  evaluate 
horizontal  bifilar  designs.  The  Sample  Applications  section 
illustrates  post-processing  of  vc  and  \  responses.  A  formula  is 
given  (see  equation  (101))  expressing  the  bifilar  mass  (N-l)  Q 
displacement  in  terms  of  yc  and  vs  to  permit  comparisons  to  be 
made  with  measured  bifilar  mass  accelerations. 

Vertical  Linear  Bifilar  BF2 

Vertical  linear  bifilar  BF2  represents  a  bifilar  vibration  ab¬ 
sorber  designed  to  reduce  N  per  rev  vibrations  perpendicular  to 
the  plane  of  rotation  of  a  helicopter  rotor.  The  bifilar  absorber 
is  assumed  to  consist  of  N  arms  equal  in  number  to  the  rotor 
blades.  Each  arm  holds  a  pendular  mass  free  to  flap  up  or  down. 
The  arms  are  assumed  to  be  equally  disposed  and  each  holds  an 
identical  absorber  mass.  The  bifilar  absorber  rotates  at  the 
helicopter  rotor  speed  fi.  Figure  9  is  a  schematic  one  arm  of  the 
bifilar.  The  derivation  of  the  equations  of  motion  employed  in 
the  BF2  model  is  similar  to  the  derivation  described  in  the 
previous  section  for  bifilar  BF1 .  Each  bifilar  mass  flapping 
displacement  is  expressed  as  the  sum  of  a  bifilar  flapping  collec¬ 
tive  and  cyclic  displacements. 

Pn  =  P0/N  +  (Pccos  ^n  +  Pssin  'i>n)  (2/N)  (77) 

Substitution  of  equation  (77)  into  the  equations  of  motion  (see 
previous  section  for  a  similar  discussion)  yields  a  nine  co¬ 
ordinate  constant  coefficient  equation  set.  The  nine  coordinates 
consist  of  bifilar  mass  collective  (PQ)  sine  cyclic  (Pg),  and 
cosine  cyclic  (Pc)  amplitudes,  and  six  hub  displacements  shown  in 
Figure  9.  Restrictions  on  the  configurations  that  can  be  analyzed 
are  bifilars  which  have  three  or  more  equally  disposed  and  identi¬ 
cal  masses.  Solution  modes  that  may  be  used  with  BF2  are  FR1, 
EG2 ,  and  TH1 .  Initial  values  of  bifilar  mass  collective,  and 
cyclic  amplitudes  may  be  specified  for  transient  analysis. 
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N  MASSES,  ONLY  ONE  MASS  SHOWN 
/3  =/3Q/N  +  (£c  COS*+£s  SIN*)(2/N) 


Figure  9.  Vertical  Linear  Bifilar  BF2 . 
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Vertical  bifilars  are  tuned  to  eliminate  hub  vibrations  induced  by 
forces  referred  to  the  rotating  system  at  a  frequency  equal  to  NV  . 
These  forces  transform  to  NQ  forces  resolved  to  non-rotating  axes. 
Post-processing  of  response  results  obtained  from  the  system 
solution  may  be  required  to  evaluate  vertical  bifilar  designs. 
The  approach  is  analogous  to  that  described  for  the  horizontal 
bifilar  (see  previous  section). 

Vibration  Isolator  IS1 

Vibration  isolator  IS1  is  a  nodal  isolator  which  may  be  placed 
between  an  excitation  source,  like  the  rotor/transmission  system, 
and  an  airframe,  to  minimize  vibrations  on  the  airframe.  The 
concept  of  the  nodal  isolator  is  described  in  Reference  7.  The 
isolator  employs  a  pivoted  bar  weight  to  balance  spring  and 
inertia  forces  to  create  a  node  or  point  of  zero  vibratory  motion 
at  the  airframe  attachment.  When  damping  is  present,  vibrations 
are  minimized  at  the  airframe  attachment  at  the  excitation  fre¬ 
quency.  In  coupled  rotor/airframe  systems  the  excitation  fre¬ 
quency  is  N  per  rev  where  N  is  the  number  of  rotor  blades.  Figure 
10  illustrates  the  basic  nodal  isolator,  whose  coordinates  are 
upper  and  lower  mass  displacements.  The  antiresonant  bar  dis¬ 
placement  is  a  function  of  upper  and  lower  mass  displacements. 

The  upper  mass  simulates  the  rotor/transmission  system.  The  lower 
mass  simulates  the  airframe.  Design  variables  are  spring  rate, 
antiresonant  bar  weight,  bar  weight  overhang,  and  distance  between 
bar  pivots.  Variables  are  selected  to  achieve  minimum  vibrations 
at  a  point  on  the  airframe  at  a  desired  excitation  frequency 
applied  by  rotor/transmission  A.  Design  variables  are  selected  to 
protect  the  isolator  from  resonance  at  the  rotor  excitation 
frequency.  Reference  7  describes  the  dynamic  attributes  of  the 
system . 
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Nodal  isolator  IS1  generalizes  the  isolator  concept  to  introduce 
six  translational  and  rotational  displacements  at  each  upper  and 
lower  isolator  attachment.  The  antiresonant  bar  weight  is 
allowed  to  translate  in  three  directions.  Figure  11  illustrates 
the  physical  isolator  and  a  schematic  corresponding  the  physical 
isolator.  Figure  11  does  not  show  a  second  antiresonant  bar 
weight  included  in  the  mathematical  model.  The  lower  bar  weight 
permits  additional  design  flexibility  and  is  shown  in  Figure  12. 
The  isolator  transmission  attachment  (node  1)  in  Figure  11  corres¬ 
ponds  to  A  in  Figure  10.  The  isolator  airframe  attachment  (node 
2)  corresponds  to  B  in  Figure  10.  Six  coordinates  at  each  attach¬ 
ment  added  to  thj.ee  coordinates  for  each  antiresonant  bar  weight 
sum  to  a  total  of  eighteen  coordinates  for  IS1.  Connection  node 
numbers  and  Euler  angles  defining  the  isolator  bar  orientation  are 
associated  with  each  attachment.  These  attributes  enable  the 
isolator  bars  to  be  connected  and  assembled  with  any  other  sub¬ 
structures  of  the  Base  Program .  A  subsequent  section  (see  Sample 
Applications)  describes  a  configuration  assembled  from  four 
isolator  bars  of  type  IS1  to  create  a  "six  degree-of-freedom" 
nodal  isolator.  Nodal  isolator  IS1  can  be  specialized  to  the 
isolator  shown  in  Figure  10. 

Rotor  Aeroelastic  Substructure  RE 3  and  External  Program  G400 

Rotor  aeroelastic  substructure  RE 3  is  used  to  represent  a  multi¬ 
blade  aeroelastic  rotor  in  forward  flight  or  hover.  The  rotor 
representation  in  this  substructure  is  defined  by  a  rotor  hub 
excitation  force  and  moment  vector  and  a  rotor  impedance  matrix. 
The  force  and  moment  vector  and  impedance  matrix  may  be  generated 
by  the  G400  aeroelastic  analysis.  The  Base  Program  combines  the 
aeroelastic  rotor  representation  with  the  impedance  elements  of 
other  dynamical  components  to  form  a  coupled  system.  The  de¬ 
finitions  of  the  rotor  hub  force  excitation  vector  and  rotor 
impedance  matrix  are  given  in  a  preceding  section  (Assembly  of 


64 


Rotor  Impedance  Matrix).  The  coordinates  represented  in  the  RE3 
substructure  are  the  six  hub  displacements  listed  in  equation 
(59).  The  RE  3  rotor  impedance  substructure  is  limited  to  the 
steady  state  forced  response  solution  mode,  FR1 .  The  frequency  of 
force  excitation  specified  in  FR1  must  be  the  same  as  that  speci¬ 
fied  for  G400  to  derive  the  rotor  impedance  matrix.  This  fre¬ 
quency  must  be  an  integer  multiple  of  rotor  speed  and  blade  number 
(w^.  =  nNfi )  .  Component  RE3  may  be  flagged  to  calculate  Higher 
Harmonic  Control  (HHC)  angles  to  minimize  vibrations,  such  as 
described  in  a  later  section  (  Higher  Harmonic  Control  for 
Vibration  Reduction).  Optimal  control  angles,  and  blade  stresses 
resulting  from  hub  displacements  and  HHC  control  inputs  may  be 
calculated.  This  is  achieved  by  utilizing  data  from  G400  defining 
the  effects  on  hub  loads  and  blade  stresses  of  HHC  angle  pertur¬ 
bations  and  hub  displacements.  Table  5  summarizes  G400 
capabilities . 

Component  RE3  may  also  be  used  to  represent  the  impedance  of  a 
substructure  which  is  not  a  rotor,  such  as  an  empennage.  Aero¬ 
dynamic  dampings  may  be  represented  in  the  impedance  matrix. 
Structural  and  inertial  contributions  to  the  impedance  matrix  may 
also  be  included. 


HIGHER  HARMONIC  CONTROL  FOR  VIBRATION  REDUCTION 


Higher  harmonic  control  (HHC)  is  a  concept  for  helicopter  vibra¬ 
tion  reduction  which  employs  high  frequency  blade  pitch  controls 
to  reduce  helicopter  vibrations.  Vibration  reduction  is  accom¬ 
plished  by  modification  of  vibratory  forces  and  moments  acting  on 
the  rotor  blades,  such  that  the  resultant  hub  forces  and  moments 
are  reduced.  This  section  describes  the  approach  used  in  the 
Coupled  Rotor/Airframe  Vibration  Analysis  to  calculate  helicopter 
vibrations  for  specified  HHC  angles  and  to  calculate  optimal  HHC 
control  angles  for  minimum  vibrations. 
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TABLE  5  -  CAPABILITIES  OF  G400  ROTOR  AEROELASTIC  ANALYSIS 


Rotor  Type 

Conventional  articulated 
Composite  bearingless 
Teetering 

Blade 

.  Redundant  torque  tube  flex  beam 
Pendulum  absorbers 
Structural  twist 
Aerodynamic  twist 

Aerodynamics 

Static  aerodynamic  data 
Unsteady  data 

Nonuniform  rotor  induced  inflow 

Stress  Calculation 
.  Mode  deflection 

External  force  integration 

Mathematical 

Uncoupled  normal  modes 
Time  historv  solution 

Hub  force  and  moment  excitation  vector  and 

rotor  impedance  matrix 

Higher  harmonic  control  matrices 
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Higher  harmonic  control  may  be  achieved  by  superimposing  non¬ 
rotating  swashplate  motions  at  the  blade  passage  frequency  upon 
the  collective  and  cyclic  control  inputs  at  0.  The  blade  passage 
frequency  is  NO  for  an  N  bladed  rotor.  The  NO  cyclic  motions  of 
the  non-rotating  swashplate  transform  to  (N-1)0  and  (N+1)0  blade 
pitch  motions  referred  to  the  rotating  system.  The  NQ  collective 
motion  of  the  swashplate  results  in  NO  blade  pitch  changes  in  the 
rotating  system.  The  harmonic  coefficients  of  the  blade  pitch 
angle  referred  to  the  rotating  system  may  be  expressed  as 


|6N-1  ■ 

’  ®N  '  0N+1}C 

(78) 

<6N-1  ■ 

’  ®N  '  0N+1} 

(79) 

The  elements  of  the  vectors  in  equations  (78)  and  (79)  are  the 
coefficients  of  cosine  and  sine  terms  in  the  pitch  angle  ex¬ 
pression  with  frequencies  (N-1)0,  NO,  and  (N+1)Q.  The  control 
inputs  induce  rotor  forces  and  moments  at  NO  in  the  non-rotating 
system. 


Higher  harmonic  control  angles  are  perturbed  in  the  G400  analysis 

Matrix  h^j  relates  the  rotor  hub 
reaction  forces  to  HHC  control  perturbations,  and  enables  HHC  to 
be  introduced  in  the  Coupled  Rotor/Airframe  Analysis.  Successive 
perturbing  of  HHC  controls  enables  columns  of  the  transfer  matrix 
to  be  identified.  The  procedure  is  similar  to  the  determination 
of  the  rotor  impedance  matrix  leading  to  equation  (60).  The 
resulting  expression  for  the  harmonics  of  hub  reaction  force  at 
frequency  NO  is 


to  obtain  a  transfer  matrix 
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(81) 
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Equation  (80)  is  the  same  as  equation  (60)  with  an  additional  term 
involving  the  HHC  angles.  Matrix  h^  has  dimensions  12x6.  Each 
submatrix  h,,„,  h, „rW  h 
6x3  . 

shown  by  equations  (78)  and  (79). 


11H'  12H'  21H' 

Each  of  the  control  vectors  6 


h,.-..  in  equation  (81)  has  dimension 


and  0  has  dimensions  3x1  as 
c  s 


Assembly  of  the  rotor  hub  excitation  forces  induced  by  HHC  with 
the  properties  of  other  substructures  is  achieved  by  utilizing 
substructure  RE3 .  The  assembly  procedure  parallels  the  hub  force 
and  moment  assembly  described  in  a  preceding  section  (see  Assembly 
of  Rotor  Impedance  Matrix).  The  resulting  system  equation  ex¬ 
pressed  in  terms  of  dependent  coordinates  is 
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Terms  in  equation  (81)  are  the  same  as  the  terms  in  equation  (64) 
except  for  the  HHC  terms.  Matrix  H11D  in  equation  (82)  is 
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Equation  (85)  is  the  same  as  equation  (27)  derived  in  a  preceding 
section  (Assembly  of  Rotor  Impedance  Matrix).  Equation  (85) 
yields  the  system  response  including  a  rotor,  without  HHC . 
Equation  (84)  yields  the  system  response  when  HHC  is  specified. 
It  is  essentially  the  same  as  equation  (1)  in  Reference  11. 
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Matrix  T  is  the  transfer  matrix  relating  HHC  inputs  to  vibratory 
response.  The  dimensions  of  the  matrices  introduced  in  equation 
(84)  are:  T  =  2Nx6 ,  A  =  2Nx2N,  H  =  2Nx6,  H^,  H  ,  H  ,  H22  =  Nx3 
where  N  is  the  number  of  coordinates  in  the  independent  coordinate 
system . 

To  derive  the  control  angle  vector  for  minimum  vibrations,  a 
scalar  performance  index  J  is  defined. 


J  = 


T 


+  {9 


c 


W„ 


(90) 


The  control  vector  must  be  chosen  to  minimize  J.  The  vector  of 
response  variables  Y,  whose  cosine  and  sine  components  are  Yc  and 
Y  ,  may  be  identified  with  the  displacement  vector  Xp,  the  inter¬ 
face  force  vector  R^,  or  a  vector  composed  of  both  displacements 
and  interface  forces.  In  this  application,  Y  is  selected  as 


Y  =  Xp  (91) 

which  identifies  Y  as  the  dependent  coordinate  vector.  The 
diagonal  weighting  matrix  Wz  permits  the  weighting  of  some 
responses  more  than  others.  Diagonal  elements  may  be  selected  to 
be  zero  except  at  a  point  on  the  structure  whose  displacements  are 
to  be  minimized.  A  typical  requirement  is  the  weighting  of 
cockpit  or  tail  responses  for  minimization.  Weighting  matrix  Wz 
has  dimensions  x  where  is  the  size  of  the  vector  of 
dependent  coordinates,  Xp. 
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The  3x3  weighting  matrix  WQ  limits  the  amplitudes  of  the  rotor 
controls.  Equation  (90)  is  analogous  to  equation  (2)  in  Reference 
LI.  The  control  state  for  a  minimum  value  of  J  is  obtained  by 
first  substituting  equation 
This  yields 

)<Vc  /  \f>Xco  / 

l<Vs\  ItXso  \ 

Equations  (91)  and  (92)  are 
express  the  performance  index  in  terms  of  control  variables.  The 
resulting  expression  is  differentiated  with  respect  to  the  control 
variables  and  equated  to  zero  to  obtain  the  control  state  for~  a 
minimum  value  of  J. 


(84)  into  equations  (67)  and  (68). 


P  J 


K  I 

)  n 


(92) 


substituted  into  equation  (90)  to 
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Equation  (93)  is  analogous  to  equation  (3)  in  Reference  11,  except 
that  it  is  expressed  in  terms  of  cosine  and  sine  components  of 
the  independent  coordinate  vector,  Xj .  An  equivalent  version  of 
equation  (93)  is  implemented  in  the  program  code  to  minimize  the 
computer  storage  requirement.  Transformations  are  applied  to  the 
T  matrix  in  the  code  as  opposed  to  the  Wz  matrix  such  as  shown  in 
equation  ( 93  )  . 
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Substituting  equation  (93)  into  equation  (92)  yields  the  minimum 
vibration  vector.  The  optimal  control  solution  and  corresponding 
vibrations  are  calculated  in  the  Base  Program.  File  data  are  read 
to  acquire  the  transfer  matrix  h^  tramsmitted  from  G400 .  A 
non-linear  dependence  of  h^  on  the  control  state  is  approximated 
by  deriving  h^  for  control  perturbations  centered  at  prior  solu¬ 
tions  obtained  from  the  Base  Program.  The  HHC  solution  is  refined 
by  iterating  between  the  Base  Program  and  G400. 

Blade  stresses  and  push  rod  loads  are  monitored  for  the  effects  of 
rotor  hub  displacement  and  higher  harmonic  control  inputs.  This 
is  achieved  by  expanding  the  harmonic  coefficient  of  a  stress  in 
terms  of  hub  displacements  and  HHC  control  angles,  such  as  was 
done  for  the  harmonics  of  hub  reaction  force  in  equation  (80). 
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[SnH^ 


*  ..  \ 


n6 


! ( 
)„  \ 


(94) 


Stress  coefficient  Sn  is  the  n-th  harmonic  of  stress  at  a  point  on 

the  rotor  blade  where  n  is  an  integer  multiple  of  rotor  speed,  0. 

Stress  coefficient  Sn  represents  either  the  flatwise,  edgewise, 

torsion  stress,  or  push  rod  load.  Stress  harmonic  coefficient  Sn 

is  a  scalar.  The  coefficient  (S  )  is  the  n-th  harmonic  coeffic- 

n  o 

ient  of  stress  without  hub  displacement  nor  HHC  controls.  Row 

matrix  [S  „]  is  a  1x12  matrix  of  coefficients  introducing  the 
rixi 

effects  of  hub  displacement.  Row  matrix  [S  ]  is  a  1x6  matrix  of 

coefficients  introducing  the  effects  of  HHC  controls.  Scalar 

(S  )  and  matrices  [S  „]  and  [S  „]  are  calculated  by  G400  and  are 
'no  nH  nO 

transmitted  to  the  Base  Program.  The  Base  Program  calculates  the 
harmonics  of  stress  or  push  rod  load  from  equation  (94).  These 
harmonics  are  used  to  calculate  the  variation  of  stress  or  push 
rod  load  with  time.  Peak  to  peak  values  are  determined  in  the 
Base  Program. 


SAMPLE  APPLICATIONS 


A  discussion  is  given  below  of  sample  applications  performed  with 
the  Base  Program  to  demonstrate  the  capabilities  of  the  analysis 
package.  The  set  of  applications  exercises  each  component 
(substructure  or  solution  mode).  The  problems  were  selected  for 
physical  significance  and  ability  to  test  program  capabilities  and 
component  combinations.  Table  6  summarizes  the  attributes  of  the 
Base  Program  sample  problems.  A  brief  description  is  keyed  to  a 
list  of  components  employed  for  the  case. 

Base  Program  Applications 

Case  1  was  used  to  verify  the  functioning  of  the  horizontal  linear 
bi filar  BF1  and  components  GF1 ,  MSI  and  FR1 .  Attributes  of  the 
UH-60  system  were  used.  The  configuration  consists  of  a  bifilar 
with  four  pendulum  masses,  coupled  to  a  modal  fuselage  with  five 
modes.  Shake  test  data  were  the  source  of  modal  properties.  A  40 
inplane  force  is  applied  to  the  hub  to  simulate  the  effects  of  a 
rotor.  The  components  of  the  applied  force  are  resolved  to  the 
fixed  system  and  are  equal  to 


fxl 

=  1700 

cos  40 t. 

N 

(lb) 

(95) 

fyi 

=  1100 

sin  40 t. 

N 

(lb) 

(96) 

The  above  values  of  the  force  coefficients  are  input  to  GF1 . 
Figure  2  shows  the  positive  directions  of  the  forces.  These 
forces  are  equivalent  to 


fx2 

=  1400 

cos  30t  + 

300 

cos  50t, 

N 

(lb) 

(97) 

fy2 

=  1400 

sin  30t  - 

300 

sin  50t, 

N 

(lb) 

(98) 
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referred  to  the  rotating  system.  The  fixed  system  4  fi  forces 

result  from  rotating  system  forces  with  frequencies  30  and  50 

referred  to  the  rotating  system.  The  rotating  system  forces 

applied  by  the  rotor  were  found  from  measured  bi filar  arm  ac¬ 
celerometer  readings.  Hub  displacements  were  assumed  to  be 

negligible  as  a  result  of  the  bifilar  having  minimized  the  hub 
motions.  The  steady  state  forced  response  of  the  system  was 

obtained  from  solution  mode  FR1 .  The  coupled  system  solution  has 
cosine  and  sine  components  (see  equation  (24))  and  the  bifilar  arm 
displacements  are  therefore 

V s  =  (Ys>c  cos  4fit  +  (Ys)s  sin  4fit,  rad  (99) 

Y c  =  (Yc)c  cos  40t  +  (YC)S  sin  40t'  rad  (100) 

Substitution  of  equations  (99)  and  (100)  in  equation  (74)  yields 

/y/3n  =  1  [[-( Ys)c)  +  (Yc)s]2  -  [<YS)S  +  (Yc)c]2)H,  rad  (101) 


For  the  bifilar  mass  30  amplitude  additional  cases  were  run  to 
find  the  bifilar  mass  amplitude  variation  with  airspeed,  shown  in 
Figure  13.  Results  are  compared  with  test  data  obtained  from  a 
bifilar  mass  mounted  accelerometer.  The  agreement  of  theoretical 
results  and  test  data  shows  the  correctness  of  the  substructure 
coupling  procedure,  the  equations  of  motion  for  the  components 
employed  in  the  configuration,  and  the  system  properties  deter¬ 
mined  from  test  data.  The  analytical  model  is  not  suitable  for 
bifilar  or  fixed  system  design  when  rotor  hub  loads  data  are  not 
available  or  at  bifilar  off-design  conditions  with  large  hub 
displacements  when  rotor  impedance  effects  are  important.  A 
suitable  general  model  for  design  studies  would  employ  an  aero- 
elastic  rotor  model,  such  as  RE3 . 
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3/REV  BIFILAR  MASS 
MOTION, G‘s 


A  FLIGHT  TEST 


Figure  13.  Bifilar  Mass  Amplitude  Variation  with  Airspeed 


Case  2  is  a  verification  of  a  vertical  linear  bifilar  BF2 .  A 
vertical  rotor  induced  force  is  applied  to  the  hub.  In  other 
respects,  the  case  is  similar  to  Case  1. 

Case  3  is  a  verification  of  beam  model  BM1 ,  constraint  component 
CN1 ,  and  eigensolution  mode  EG1.  The  case  yields  the  mode  shapes 
and  frequencies  for  a  five  element  uniform  beam.  The  beam  is 
constrained  from  motion  in  all  directions  except  vertical  trans¬ 
lation.  Figure  14  shows  mode  shapes  and  frequencies  for  a  ten 
element  beam  model  similar  to  Case  3.  Results  are  in  close 
agreement  with  classical  results.  The  ten  element  beam  case 
required  approximately  15K1Q  decimal  working  storage  words  (1K10  = 
1024  words),  which  together  with  the  storage  required  for  instruc¬ 
tions  and  common  blocks  (17K1Q  decimal  words)  shows  that  this 
problem  has  approached  the  32K^0  decimal  word  limit  for  inter¬ 
active  CDC  system  usage  (see  Reference  17).  The  beam  model  BM1 
could  have  been  used  with  any  solution  mode.  In  particular,  a 
forced  response  study  with  FR1  and  rotor  system  RE 3  could  have 
been  performed  to  study  the  vibration  behavior  of  a  helicopter 
preliminary  design. 

Case  4  is  a  verification  of  the  fixed  absorber  FA1 ,  with  parameter 
variation  mode  PV1 .  A  mass  is  attached  to  the  base  of  the  ab¬ 
sorber  by  means  of  an  MSI  component  to  simulate  a  two  mass  system 
joined  by  a  spring  and  damper.  An  oscillating  force  is  applied  to 
the  base  of  the  absorber  and  a  steady  state  forced  response 
solution  is  obtained  The  forcing  frequency  is  specified  to  range 
between  one  and  ten  Hz.  This  parameter  variation  is  achieved  by 
means  of  input  to  PV1 .  Vibratory  response  results  agree  with 
analytical  derivations.  Analytical  results  were  derived  manually 
by  application  of  the  Laplace  transform  method  to  the  assumed 
system . 
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Case  5  is  similar  to  Case  4  except  that  a  time-history  solution 
mode  TH1  substitutes  for  steady-state  forced  response  mode  FR1  and 
parameter  variation  mode  PV1 .  The  case  is  also  used  to  verify  the 
general  (GEN)  control  mode  line  printer  output  suppression  capa¬ 
bility.  The  case  transmits  data  to  a  plot  file  and  the  Plot 
Package  is  used  to  plot  results.  Figure  15  illustrates  the 
variation  with  time  of  the  responses  of  the  two  masses.  The 
steady-state  values  of  vibratory  response  agree  with  the  values 
obtained  from  Case  4,  and  with  analytical  derivation.  A  rigid 
body  component  of  displacement  of  the  two  mass  system  shows  a 
linear  increase  with  time  and  agrees  with  analytical  derivation. 

Case  6  is  the  same  as  Case  5  except  that  it  is  used  to  check  the 
restart  option  in  the  time-history  solution  mode.  The  case 
employs  time*.  ! story  restart  data  transmitted  to  a  file  from  a 
prior  time  history  run. 

Case  7  consists  of  one  nodal  isolator  of  type  IS1  coupled  to  upper 
and  lower  masses,  connected  by  a  spring  such  as  shown  in  Figure 
10.  The  isolator  is  free  to  move  vertically.  The  frequency 
response  of  the  system  co  a  force  applied  to  the  isolator  upper 
node  is  shown  in  Figure  16.  The  isolator  capability  to  create  a 
point  of  zero  vibration  at  the  lower  attachment  node  is  evident  at 
the  antiresonant  frequency  (27  Hz). 

Cases  8  to  10  serve  as  a  systematic  synthesis  and  verification  of 
features  employed  in  a  hypothetical  six  degree-of-freedom 
isolator  model,  Case  11. 

Case  8  is  similar  to  Case  7  in  the  character  except  that  the 
isolator  is  free  to  move  vertically  and  laterally.  Two  modal 
structures  (MSI)  were  connected  to  the  isolator  at  nodes  one  and 
two.  A  force  of  type  G FI  was  applied  to  node  1  of  the  isolator. 
The  configuration  is  shown  in  Figure  17. 
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Figure  15.  Variation  with  Time  of  the  Dispiacment  of  a  Two-mass 
System  to  which  an  Oscillating  Force  is  Applied. 
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Figure  16.  Frequency  Response  of  a  Simple  Nodal  Isolator 
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Figure  17.  Configurations  Used  to  Test  I  SI  Component. 
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Case  9  is  similar  to  case  8  except  that  the  GF1  force  was  applied 
to  node  3  of  an  MSI  substructure  (see  Figure  17)  which  connects 
the  force  application  point  and  the  substructure  attachment  to  the 
upper  isolator  node  (node  2).  The  upper  MSI  substructure  simu¬ 
lates  a  transmission  placed  between  a  rotor  and  an  isolator.  The 
lower  MSI  substructure  simulates  a  fuselage. 

Case  10  generalizes  Case  9  to  a  configuration  employing  six  modes 
for  transmission  representation,  six  modes  for  fuselage  representa¬ 
tion,  and  an  isolator  of  type  IS1.  All  the  modes  are  rigid-body 
modes.  A  steady  state  forced  response  solution  is  obtained. 

Case  11  consists  of  a  six  degree-of-freedom  isolator  constructed 
from  four  isolator  bars  of  type  IS1  connected  at  right  angles  to 
each  other.  Six  rigid-body  transmission  modes,  and  six  rigid-body 
fuselage  modes  are  simulated.  Six  hul  forces  and  moments  are 
applied  to  the  upper  transmission  to  simulate  the  loading  applied 
by  a  rotor.  Isolator  characteristics  were  derived  from  a  half 
scale  UH-60  isolator.  An  illustration  of  the  model  is  shown  in 
Figure  18.  A  total  of  nine  connection  nodes  is  used.  This  case 
introduces  substructure  Euler  angles  to  check  the  correct  forma¬ 
tion  of  angular  transformations  from  local  to  global  axes.  Bar 
Euler  angles  of  0,  90,  180  and  270  degrees  are  introduced.  Figure 
19  shows  the  frequency  response  of  the  isolator  at  a  fuselage 
attachment.  The  variations  with  forcing  frequency  of  three 
components  of  translation  are  shown  at  this  point.  Figure  20  is  a 
three  dimensional  view  of  vertical  displacements  at  the  four 
fuselage  attachments  (see  Reference  18  for  plot  types  that  can  be 
requested ) . 

Case  12  verifies  the  ability  of  the  program  to  employ  rotor 
aeroelastic  model  RE2  (E927  rotor).  A  five  mode  fixed  system  was 
coupled  to  the  rotor.  The  rotor  employs  twelve  coordinates 
consisting  of  flap  collective,  lead  collective,  two  flap  cyclic 
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Figure  18.  six  Degree-of-Freedom  Isolator 
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Figure  19. 


Frequency  Response  of  Six  Degree-of- 
Isolator  at  a  Fuselage  Attachment. 


90 


Figure  20.  Six  Degree-of-Freedom  Isolator  Vertical 
Responses  at  Fuselage  Attachment  Points. 


modes,  two  lead  cyclic  modes,  two  bending  collective  modes,  and 
four  bending  cyclic  modes.  An  eigensolution  of  the  system  was  ob¬ 
tained  using  the  eigensolution  mode  with  damping,  EG2 .  Results 
agreed  with  results  obtained  from  the  unmodified  E927  analysis, 
which  included  an  elastic  fuselage. 

Case  13  verifies  the  ability  of  the  program  to  employ  rotor 
aeroelastic  analysis  RE3  and  to  serve  as  a  tool  to  study  effici¬ 
ently  the  effects  of  variations  in  components  internal  to  the  fuse 
Program.  A  rotor  hub  load  excitation  vector  and  a  hub  impedance 
matrix  were  transmitted  by  G400.  A  hypothetical  full  scale  four 
blade  rotor  and  four  fixed  system  modes  were  specified.  One  flap, 
one  lead,  three  uncoupled  flatwise,  two  uncoupled  edgewise,  and 
one  uncoupled  torsion  blade  modes  were  used.  Airspeed  of  120 
knots  was  specified,  and  uniform  inflow  was  used.  Figures  21  and 
22  show  the  effects  of  changes  in  the  frequency  of  the  third 
fuselage  mode  on  hub  loads  and  accelerations ,  respectively  (magni¬ 
tudes  of  leads  and  accelerations  should  not  be  interpreted  as 
being  representative  of  an  actual  helicopter).  Selected  elements 
of  the  rotor  impedance  matrix  were  derived  analytically  and  agreed 
with  the  G400  results.  Verification  of  the  ability  of  program 
package  to  derive  and  use  rotor  hub  loads  and  impedances  for 
rotors  employing  variable  inflow  is  covered  in  the  subsequent 
text.  Case  13  verifies  the  END  CASE  command  (see  Reference  17). 
The  first  set  of  data  in  Case  13  employs  a  flag  in  RE 3  to  display 
the  output  solution  vector  in  feet  and  the  rotor  hub  interface 
forces  (see  equation  (60)).  The  data  following  the  END  CASE 
command  request  the  output  to  be  displayed  in  units  of  accelera¬ 
tion  of  gravity  (see  equations  (28)  and  (29)). 

EXTERNAL  PROGRAMS  IN  THE  PROGRAM  PACKAGE 

Descriptions  are  given  here  of  the  functions  and  assumptions  of 
the  External  Programs  available  in  the  Coupled  Rotor/Airframe 
Vibration  Analysis  package . 
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Figure  22.  Effects  on  Rotor  Hub  Accelerations  of  Variations 
in  Third  Fixed  System  Mode  Frequency. 
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Rotor  Aeroelastic  Analysis  E927 


The  E927  Rotor  Aeroelastic  Stability  Analysis  (see  Reference  6)  is 
a  linear  analysis  yielding  dynamic,  and  aerodynamic  mass,  damping, 
and  stiffness  matrices  of  the  rotor  system  in  hover.  These 
matrices  are  associated  with  perturbations  of  the  coordinates  of 
the  rotor  system.  Coordinates  comprise  blade  generalized  dis¬ 
placements,  both  physical  (e.g.  flap  angle)  and  modal.  Five  hub 
displacements  are  represented;  lateral,  longitudinal,  and  vertical 
translations,  and  hub  pitch  and  roll.  Linearization  occurs  about 
a  rotor  trim  state.  The  trim  vector  is  obtained  from  G400  or  any 
other  suitable  analysis.  Mass,  damping,  and  stiffness  matrices 
are  transmitted  to  a  file  and  read  by  the  Base  Program.  These  are 
assembled  with  matrices  for  other  components  selected.  The  Base 
Program  eigensolution  mode  is  used  to  evaluate  the  stability  of 
the  rotor  system  coupled  to  other  components.  Steady-state  forced 
response  may  be  calculated  if  a  rotor  hub  force  and  moment  excita¬ 
tion  vector  is  input  to  the  Base  Program.  Hub  force  and  moment 
excitation  data  are  obtained  from  the  G400  analysis  or  any  other 
suitable  analysis. 


Rotor  Aeroelastic  Analysis  G400 

The  G400  Rotor  Aeroelastic  Analysis  yields  the  time-history 
response  of  conventional  articulated,  teetering,  and  composite 
bearingless  -otors.  Options  permit  the  use  of  variable  rotor 
induced  inflow,  static  blade  aerodynamic  data,  and  unsteady 
aerodynamic  data.  Blade  mounted  pendulum  absorbers  are  modeled. 
Forward  flight  and  hover  are  represented.  The  program  calculates 
the  hub  excitation  vector  for  a  fixed  hub,  and  the  rotor  impedance 
matrix  relating  six  hub  forces  and  moments  to  harmonic  perturba¬ 
tions  of  the  hub.  These  data  are  transmitted  to  a  file  which  is 
read  by  the  Base  Program.  The  rotor  hub  excitation  vector  and  the 
hub  impedance  are  assembled  with  mass,  damping,  and  stiffness 
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elements  of  internal  components  in  the  Base  Program.  The  steady 
state  forced  response  mode  in  the  Base  Program  is  invoked  to 
obtain  coupled  system  responses. 

Rotor  aeroelastic  analysis  G400  was  created  to  overcome  inadequa¬ 
cies  of  applications  of  conventional  analyses  to  composite  bearing¬ 
less  rotors.  The  analysis  is  based  on  a  normal  modes  model  of  a 
non-linear  rotor  blade.  References  8  and  9  describe  the  basis  of 
the  method.  A  coordinate  transformation  is  applied  to  second 
derivatives  of  blade  bending  deflections  as  opposed  to  bending 
deflections.  This  transformation  enables  uncoupled  normal  modes 
to  be  used  to  achieve  non-linear  twist  representations  which 
satisfy  the  blade  root  and  tip  boundary  conditions.  Uncoupled 

modes  derive  from  eigensolutions  which  employ  uncoupled  flatwise 
bending,  edgewise  bending,  and  torsion  equations.  Up  to  four 
flatwise,  three  edgewise,  and  three  torsional  elastic  modes  may  be 
used.  Blade  articulation  reduces  the  number  of  allowed  modes.  A 
Galerkin  method  yields  modal  differential  equations  which  are 
integrated  to  obtain  the  time  history  response  of  the  blade. 
Blade  stress  calculation  options  are  mode  deflection  and  external 
force  integration  methods  (see  Reference  10).  Single  blade  loads 
are  analyzed  harmonically  to  obtain  harmonics  of  rotor  loads  for  a 
fixed  hub.  The  hub  is  perturbed  by  specified  harmonic  motions  and 
the  resulting  hub  load  perturbations  are  analyzed  harmonically  to 
obtain  the  rotor  impedance  matrix  (see  Assembly  of  Rotor  Impedance 
Matrix).  Control  perturbations  yield  matrices  required  with  the 
HHC  solution  mode  (see  Higher  Harmonic  Control  for  Vibration 
Reduction).  Table  5  lists  the  capabilities  of  the  G400  analysis. 
All  rotor  types  are  based  on  a  single  blade  analysis,  except  for 
teetering  rotors  which  are  treated  as  two  bladed  rotors.  Variable 
rotor  induced  inflow  is  introduced  by  executing  sequentially  G400 
and  the  F389  Prescribed  Wake  _>.otor  Inflow  Analysis.  Refinement 
with  variable  inflow  ~'ced  calculation  of  the  hub  excitation 
vector  and  impedance  n.c-erix.  The  same  interfacing  of  the  Base 
Program  and  G400  occurs  as  in  the  uniform  inflow  case.  Fuselage 
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induced  inflow  at  the  rotor  disc  is  introduced  by  transferring 
harmonics  of  these  data  from  WABAT  to  inflow  program  F389. 

Prescribed  Wake  Rotor  Inflow  and  Flow  Field  Analysis  F389 

Prescribed  wake  rotor  inflow  and  flow  field  analysis  F389  cal¬ 
culates  the  aerodynamic  flow  field  induced  by  a  rotor.  The 
analysis  may  be  used  to  calculate  rotor  induced  velocities  in  the 
plane  of  the  rotor  or  at  user  selected  points  -  e.g.,  at  the 
empennage . 

The  F389  analysis  is  based  on  a  lifting  line  rotor  blade  re¬ 
presentation.  A  vortex  wake  of  prescribed  shape  is  assumed  to 
trail  from  the  rotor  blades.  The  assumption  built  into  the 
program  is  that  the  wake  form  is  a  skewed  helix  (classical  wake) 
as  shown  in  Figure  23.  The  skew  angle  and  distance  between 
successive  helixes  is  calculated  by  assuming  the  wake  to  be 
convected  from  the  rotor  by  a  uniform  self-induced  rotor  inflow 
derived  from  momentum  balance  considerations,  and  from  contribu¬ 
tions  from  components  of  vehicle  flight  speed.  The  distance  from 
the  blade  at  which  the  trailing  vortices  roll  up  into  a  single  tip 
vortex  may  be  specified.  The  wake  is  terminated  at  a  finite 
number  of  wake  revolutions.  The  method  of  finding  unknown  circu¬ 
lations  on  blades  and  wake  is  described  in  Reference  12.  The 
state  of  motion  of  the  blades  is  assumed  to  be  known.  This  state 
is  derived  from  any  suitable  analysis  external  to  F389. 

The  procedure  consists  of  using  the  Kutta-Joukowski  condition  to 
derive  a  formula  expressing  blade  circulation  as  the  product  of 
blade  section  lift  curve  slope,  and  blade  section  effective  angle 
of  attack.  The  effective  angle  is  the  geometric  angle  of  attack 
corrected  for  inflow  induced  by  rotor  blades  and  wake.  Circula¬ 
tion  continuity  (Helmholtz  law)  is  used  to  express  the  circulation 
values  in  the  trailing  wake  in  terms  if  the  circulations  on  the 
blades.  A  steady-state  of  flight  is  assumed,  and  as  a  result,  the 
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Figure  23.  Computer  Wake  Representation-Undistorted  Wake 


azimuthal  variation  of  circulation  is  periodic.  Flow  velocities 
induced  by  blades  and  trailing  vortex  wake  are  found  from  the 
Biot-Savart  law  and  contain  circulation  variables  as  unknowns. 
The  induced  velocity  expressions  are  substituted  in  the  formula 
for  blade  circulation  to  obtain  a  finite  set  of  linear  simultan¬ 
eous  equations  which  are  solved  for  the  unknown  circulations. 

After  the  circulations  are  known,  induced  velocities  can  be 
calculated  anywhere  in  the  flow  field  by  application  of  the  Biot 
Savart  law. 

Section  airfoil  lift  coefficients  are  approximated  by  straight 
line  curve  fits  to  airfoil  data,  approximating  unstalled  and 
stalled  regions  of  lift  curves  (see  Figure  24). 

The  F389  program  can  be  made  to  calculate  the  flow  field  consis¬ 
tent  with  the  motions  of  an  elastic  blade.  This  is  achieved  by 
calculating  the  component  of  airflow  normal  to  the  median  surface 
of  the  rotor  blade  induced  by  the  motion  of  the  elastic  blade, 
which  excludes  the  rotor  induced  flow  contribution  to  the  airflow. 
This  airflow  component  is 

w 1  =  UT0  +  Up  -  ROTOR  INDUCED  FLOW  (102) 

The  calculation  of  w.^  is  performed  in  G400.  The  angle,  0,  is  the 
blade  pitch  angle.  The  variable,  w^,  is  automatically  transmitted 
to  F389  by  G400.  Program  F389  adds  the  rotor  induced  inflow  to  w^ 
to  calculate  the  effective  angle  of  attack  of  a  blade  section. 
Program  F389  can  represent  an  elastic  rotor  since  no  assumptions 
are  made  restricting  UT  and  to  non-elastic  blades.  It  should 
be  noted  that  anoth'  mode  of  F389  operation  is  to  employ  the 
flapping  harmonics  to  define  the  blade  motion  input.  This  mode  of 
calculation  is  less  accurate  than  the  mode  employing  w^  because  it 
omits  elastic  displacement  contributions  to  airflow  velocity 
components . 
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Figure  25.  Calculation  Sequences  For  The  Coupled 
G400/F389  Programs 
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Determination  of  rotor  induced  flow  requires  a  coupling  of  F389 
and  a  rotor  aeroelastic  analysis.  In  the  vibration  analysis 
package,  this  is  achieved  by  transmitting  data  automatically 
between  F389  and  G400.  Data  transmission  and  the  invoking  of  F389 
and  G400  is  controlled  by  a  Fortran  driver  program,  VIBIT.  Driver- 
program  VIBIT  causes  G400  to  be  run  first.  Uniform  rotor  induced 
inflow  is  used.  Program  G400  transmits  w^  to  F389  which  calcu¬ 
lates  the  rotor  induced  flow.  The  induced  flow  is  transmitted  by 
F389  to  G400  which  recalculates  the  blade  aeroelastic  response. 
The  execution  sequence  G400/F389/G400  is  referred  to  as  one  cycle 
of  operation.  An  arbitrary  number  of  cycles  may  be  performed 
until  the  user  is  satisfied  that  consistency  of  blade  circulations 
(lifts),  blade  motions,  and  aerodynamic  velocities  is  achieved 
among  the  programs.  Figure  25  illustrates  the  V1BIT/G400/F389 
calculation  sequences.  After  a  converged  solution  is  achieved 
calculations  may  be  performed  in  the  same  run  to  obtain  rotor 
output  data  -  e.g.,  hub  excitations  and  impedance. 

The  F389  program  may  be  used  to  calculate  rotor  induced  velocities 
at  the  empennage  to  provide  data  for  the  calculation  of  rotor  wake 
induced  empennage  excitations  by  the  RIEVA  program  (see  Figure  1). 
The  empennage  flow  field  is  assumed  to  have  no  effect  on  the  rotor 
and  no  iterative  calculations  involving  F389  and  RIEVA  are  per¬ 
formed  . 

The  F389  program  accepts  input  defining  velocities  induced  by  a 
fuselage  or  gust.  The  velocity  field  is  defined  as  harmonics  of 
inflow  at  the  rotor.  Fuselage  flow  field  data  may  derive  from 
WABAT  (see  Figure  1).  No  change  in  fuselage  velocities  is  assumed 
to  occur  as  a  result  of  the  rotor  wake  influence  on  the  fuselage. 

Wing  and  Body  Aerodynamic  Technique  (WABAT) 

The  Wing  and  Body  Aerodynamic  Technique  (WABAT)  is  used  to  calcu¬ 
late  a  fuselage  flow  field.  The  flow  velocities  may  be  obtained 
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on  the  surface  of  the  fuselage  or  at  arbitrary  (off-body)  points. 
The  velocity  components  in  the  plane  of  the  rotor  can  be  expressed 
as  a  harmonic  series  of  blade  azimuth  angle  whose  coefficients 
depend  on  blade  radius.  The  harmonic  coefficients  are  transmitted 
to  F389  to  define  the  effects  of  fuselage  induced  flow  at  the 
rotor . 

To  obtain  the  solution  to  fuselage  flow,  the  surface  of  the 
fuselage  is  divided  into  planar  quadrilateral  panels  on  which 
aerodynamic  sources  of  constant  strength  are  placed.  The  flow  is 
assumed  to  be  incompressible  -  i.e.,  the  airspeed  of  the  heli¬ 
copter  is  assumed  to  be  low.  The  sources  satisfy  the  Laplace 
equation,  which  governs  incompressible,  inviscid  (potential)  flow. 
Strengths  of  the  panel  sources  are  determined  from  the  condition 
that  the  flow  normal  to  the  fuselage  is  zero  or  that  it  is  some 
specified  value  like  that  induced  by  an  engine  intake  or  exhaust. 
The  number  of  solution  points  is  specified  to  be  equal  to  the 
number  of  sources  of  unknown  strength.  A  set  of  linear  simultan¬ 
eous  equations  for  the  panel  source  strengths  is  obtained  and 
solved.  When  the  source  strengths  are  found,  the  velocity  field 
induced  by  the  fuselage  at  any  field  point  may  be  calculated  by 
summing  the  contributions  of  the  panels.  Appendix  B  of  Reference 
13  gives  further  details  of  the  method. 

Lifting  wings  are  represented  by  distributing  horseshoe  vortices 
(see  Figure  26)  in  the  mean  surfaces  of  the  wings  and  in  the  wing 
wake  (trailing  vortices).  The  solution  technique  for  wing  and 
body  configurations  consists  of  solving  for  the  distribution  of 
lift  on  the  wing  without  the  body.  This  lifting  solution  is 
achieved  by  means  of  a  modified  Multhropp  analysis  (see  Appendix  B 
Reference  13).  The  Multhropp  lift  distribution  assigns  circu¬ 
lation  weighting  factors  among  the  horseshoe  vortices.  In  other 
words,  the  relative  strengths  of  the  circulations  among  the 
horseshoe  vortices  distributed  on  a  chord  at  each  spanwise  station 
are  known.  The  magnitudes  of  the  circulations  are  not  known.  As 
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a  result  of  assignment  of  circulation  weighting  factors  and  a 
summing  of  the  circulations  on  each  chord,  the  circulation  of  each 
vortex  can  be  expressed  as  some  fraction  of  the  total  circulation 
at  a  spanwise  station.  The  number  of  unknowns  associated  with  the 
lifting  wing  at  each  spanwise  station  is  therefore  reduced  to  one, 
which  is  the  total  circulation  at  the  station.  The  total  number 
of  lifting  wing  unknowns  is  N,  where  N  is  the  number  of  spanwise 
stations  on  the  wing.  The  determination  of  the  circulation 
strength  magnitudes  is  preceded  by  placing  source  panels  on  the 
wing  to  represent  the  thickness  effects  of  the  wing  on  the  flow. 
The  unknown  circulations  and  source  strengths  are  then  obtained  by 
summing  the  induced  velocity  effects  of  body  and  wing  source 
panels  and  wing  horseshoe  vortices,  and  equating  these  velocities 
to  specified  velocities  on  the  wing  and  body  surfaces.  When  a 
wing  is  present,  solution  points  are  defined  at  each  panel  point 
on  wing  and  body.  In  addition  there  are  N  solution  points  placed 
at  the  3/4  chord  of  the  aft-most  horseshoe  vortex  at  each  spanwise 
station.  The  N  additional  points  are  introduced  to  increase  the 
number  of  boundary  conditions  to  accommodate  .he  N  additional 
unknowns  (circulations)  of  the  lifting  wing.  As  a  result,  the 
number  of  unknowns  and  solution  points  are  equal  and  the  equations 
are  solved  for  source  and  horseshoe  vortex  strengths.  A  wing  and 
body  solution  satisfying  the  boundary  conditions  of  the  fuselage 
and  wing  is  consistently  obtained. 

The  synthesis  of  the  fuselage  geometry  is  facilitated  by  an 
automated  technique  of  cross-section  development  from  input  data. 
Input  data  define  a  limited  number  of  points  on  the  cross  section. 
This  is  performed  by  the  Automated  Paneling  Technique  (APT/Plot) 
program.  Appendix  A,  Reference  13  describes  the  basis  of  the 
method.  Incorporated  in  this  program  is  a  Tektronix  plot  capa¬ 
bility  which  permits  the  user  to  verify  the  correctness  of  the 
panel  input  and  transmits  the  panel  data  by  means  of  a  file  to  the 
WABAT  potential  flow  program.  The  user  may  develop  the  fuselage 
geometry  in  subsections.  Each  section  consists  of  a  number  of 
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panels.  The  Tektronix  plot  package  assembles  the  sections  into  a 
complete  fuselage  and  displays  the  resulting  fuselage.  The 
division  of  the  fuselage  into  sections  facilitates  the  development 
of  the  complete  fuselage  geometry.  Parts  of  the  fuselage  -  e.g., 
engine  nacelles  -  requiring  complex  panel  geometries  may  be 
developed  and  verified  separately.  The  plot  package  is  able  to 
provide  views  of  the  fuselage  from  any  angle  and  distance  to 
permit  the  user  to  study  and  understand  the  data.  Sections  of  the 
fuselage  may  be  separated  from  each  other  to  understand  the  data 
and  character  of  the  sections.  Figures  27  to  29  illustrate  views 
of  the  fuselage  derived  from  the  APT/Plot  program.  Reference  14 
provides  user  instructions  for  the  APT/Plot  program. 

Panel  geometry  data  are  transmitted  by  means  of  a  file  to  the 
WABAT  potential  flow  program  (Figure  1)  to  calculate  the  wing  and 
fuselage  induced  flow.  The  matrix  of  aerodynamic  influence 
coefficients  of  the  wing  and  body  singularities  (sources  and 
vortices)  may  be  saved  on  a  file  for  subsequent  applications. 
Employment  of  a  stored  influence  coefficient  matrix  reduces  the 
run  time  for  subsequent  cases  using  the  same  wing/fuselage  geo¬ 
metry  when  attitude  ( angle-of-attack,  sideslip)  is  varied.  User 
instructions  are  provided  in  Reference  14.  A  secondary  mode  of 
the  APT/Plot  package  may  be  invoked  to  obtain  a  fuselage  surface 
flow  display  after  the  WABAT  program  has  been  run.  The  arrows  in 
Figures  27  to  29  illustrate  the  flow  directions  of  perturbational 
flows  on  fuselage  panels. 

Rotor  Wake  Induced  Empennage  Vibrations  (RIEVA) 

The  Rotor  Induced  Empennage  Vibration  Airloads  (RIEVA)  program 
calculates  the  harmonics  of  empennage  lift,  rolling  moment  and 
pitching  moment  induced  by  a  rotor  wake  at  the  blade  passage 
frequency  (NQ).  These  harmonic  coefficients  are  input  to  the 
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Coupled  Rotor/Airframe  Vibration  Analysis  Base  Program  to  calcu¬ 
late  the  vibratory  response  of  a  helicopter  to  rotor  wake  excita¬ 
tions  of  the  empennage.  The  steady  state  forced  response  mode  in 
the  Base  Program  is  used  to  obtain  a  solution  for  the  coupled 
system . 

The  empennage  is  represented  as  a  three  dimensional  lifting 
surface  of  arbitrary  shape.  The  empennage  may  include  horizontal 
and  vertical  stabilizers.  The  flow  is  assumed  to  be  incompres¬ 
sible  and  unsteady.  The  wing  and  its  wake  are  represented  by  a 
doublet  lattice  on  the  mean  surface  of  the  wing  and  a  vortex 
lattice  on  the  wake.  Figure  30  shows  the  lattice  geometry.  The 
perturbation  velocities  induced  by  the  doublet  and  vortex  lattices 
are  calculated  by  the  Biot-Savart  law.  At  each  instant  of  time 
the  strengths  of  doublets  and  vortices  must  satisfy  the  boundary 
conditions  of  no  flow  normal  to  the  wing  surface,  the  Kutta 
condition  of  smooth  flow  from  the  trailing  edge,  and  continuity 
relationships  among  the  wake  vortices  and  at  the  trailing  edge 
interface  between  doublets  and  vortices.  The  distribution  of 
circulations  among  the  trailing  vortices  is  known  because  vortices 
of  known  strengths  are  convected  from  the  trailing  edge  at  the 
speed  of  the  flow.  At  the  trailing  edge  there  is  an  equivalence 
between  the  lifting  surface  doublet  strengths  and  the  strengths  of 
the  shed  vortices.  As  a  result  of  these  relationships  the  number 
of  unknowns  reduces  to  the  number  of  constant  doublet  panels  on 
the  empennage.  These  unknowns  are  assumed  to  satisfy  as  many 
simultaneous  equations  for  normal  velocities  as  there  are  boundary 
condition  points  specified  on  the  empennage.  The  doublet 
strengths  are  obtained  by  simul taneously  satisfying  these  equa¬ 
tions  at  each  time  instant.  The  wake  is  assumed  to  begin  at  the 
calculation  starting  time.  When  doublet  strengths  are  known 
empennage  surface  pressures  may  be  calculated  at  each  instant  of 
time.  The  passage  over  the  empennage  of  rotor  vortices  causes  the 
flow  over  the  empennage  to  be  unsteady.  When  the  empennage 
aerodynamic  response  is  determined  to  be  periodic  at  the  blade 
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Figure  30.  Empennage  and  Wake  Lattice  Geometry  in  RIEVA  Analysis 


passage  frequency  (NO)  the  time-stepping  calculation  is  termin¬ 
ated.  Surface  pressures  are  integrated  to  obtain  empennage  lift, 
pitching  moment  and  rolling  moments.  These  are  harmonically 
analyzed  at  frequency  NQ .  Vibratory  response  of  a  helicopter  to 
rotor  wake  induced  empennage  excitations  is  obtained  by  input  of 
these  harmonics  to  the  Base  Program.  Empennage  loads  are  assumed 
to  be  applied  to  this  fuselage  junction. 

Close  passage  of  rotor  induced  vortices  to  the  empennage  is 
approximated  by  assuming  surface  pressures  to  be  induced  by  Ran 
kine  viscous  vortices  (see  Reference  15).  These  consist  of 
vortices  with  a  solid  core  and  a  potential  vortex  region  outside 
the  core.  The  vortex  is  assumed  to  be  more  than  one  vortex  core 
radius  away  from  the  lifting  surface.  If  the  vortex  is  closer 
than  one  core  radius  to  the  lifting  surface,  it  is  moved  to  a  new 
position  such  as  described  in  Reference  15. 

The  rotor  induced  wake  is  divided  into  a  near  wake  region  and  a 
far  wake  region  such  as  shown  ir.  Figure  31.  The  near  wake  region 
is  contained  by  two  vertical  streamwise  planes  parallel  to  the 
flight  direction,  which  also  contain  the  empennage.  The  far  wake 
is  outside  this  region.  The  RIEVA  program  provides  an  internal 
calculation  of  the  near  wake  from  data  defining  the  flight  speed, 
tip  path  plane  angle,  position  of  the  empennage,  and  rotor  momen¬ 
tum  induced  flow.  The  far  wake  effects  are  represented  as  rotor 
induced  flow  velocity  harmonics.  These  harmonics  are  obtained 
from  F389  and  are  input  into  RIEVA.  The  user  may  elect  not  to 
input  far  wake  harmonics  if  these  effects  are  small.  (Omission  of 
far  wake  data  simplifies  RIEVA  program  usage,  and  has  been  found 
to  be  justified  on  the  basis  of  applications  made  to  this  date.) 

The  following  approximations  are  made  in  the  present  version  of 
the  RIEVA  analysis.  It  is  assumed  first  that  the  air  pressures 
acting  on  the  empennage  are  not  influenced  by  the  motion  of  the 
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Figure  31.  An  Example  of  Splitting  of  Rotor 


Wake  into  Two  Parts. 


empennage.  A  second  approximation  is  that  the  empennage  dis¬ 
placement  consists  of  plunge,  pitch,  and  roll  -  i.e.,  the  empen¬ 
nage  is  assumed  to  displace  as  a  rigid  body.  The  corresponding 
generalized  forces  acting  on  the  empennage  are  lift,  pitching 
moment,  and  rolling  moment.  These  are  assumed  to  act  on  the 
elastic  structure  at  the  empennage  fuselage  junction.  When  a 
modal  representation  of  the  fuselage  is  used  it  is  necessary  to 
define  the  components  of  modal  displacement  at  this  junction  - 
e.g.,  plunge,  pitch,  and  rolling  elements  of  the  modal  matrix  -  to 
calculate  the  generalized  forces  induced  by  the  empennage. 
Refinement  of  the  calculation  of  generalized  forces  can  be  accomp¬ 
lished  by  multiplying  empennage  elastic  mode  shapes  by  the  dis¬ 
played  values  of  pressures  on  the  lifting  surfaces  to  obtain 
generalized  forces.  The  code  does  not  provide  this  calculation. 

The  assumption  that  air  pressures  on  the  empennage  are  not  in¬ 
fluenced  by  motion  of  the  empennage  is  equivalent  to  neglecting 
aerodynamic  damping  in  plunge,  pitch,  and  roll,  and  aerodynamic 
stiffness  in  pitch.  These  effects  may  be  introduced  as  follows. 
Aerodynamically  induced  impedance  elements  can  be  calculated  by 
employing  an  unsteady  aerodynamic  analysis.  A  simple  approach 
would  be  to  use  Theodorsen  unsteady  aerodynamics  (see  Reference 
16)  and  a  strip  theory  assumption.  Introduction  of  the  empennage 
impedance  matrix  may  be  achieved  by  connecting  a  substructure  of 
type  RE 3  to  the  point  of  application  of  empennage  forces  (see  the 
subsection  Rotor  Aeroelastic  Substructure  RE3  and  External  Program 
G400  for  use  of  RE3  as  a  general  impedance  substructure). 

There  are  several  alternative  ways  to  representing  the  structural 
and  inertial  properties  of  the  empennage.  One  method  is  to 
provide  a  modal  representation.  Fuselage  modes  would  be  deter¬ 
mined  with  the  empennage  assumed  to  be  present.  An  alternative 
approach  is  to  calculate  fuselage  modes  without  the  empennage. 
Empennage  properties  would  be  represented  as  contributions  to  the 
impedance  matrix  of  substructure  RE3 ,  which  are  connected  to  the 
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fuselage  at  the  empennage/fuselage  junction.  Impedance  elements 
would  also  include  aerodynamic  contributions. 

A  third  approximation  in  the  analysis  is  that  it  does  not  include 
the  contribution  to  the  streamwise  velocity  component  induced  by 
vorticity  shed  from  the  rotor.  This  increment  may  be  important 
when  a  rotor  vortex  passes  near  the  vertical  stabilizer.  Com¬ 
ponents  of  rotor  induced  velocities  normal  to  the  lifting  surfaces 
and  perturbational  streamwise  components  induced  by  the  lifting 
surfaces  are  represented  and  are  essential  for  obtaining  the  RIEVA 
solution . 


Base  Program  Plot  Package 

The  Base  Program  stores  solution  information  on  file  for  pro¬ 
cessing  by  the  Base  Program  Plot  Package.  Line  printer  output 
produced  by  the  Base  Program  may  be  suppressed  optionally.  The 
Plot  Package  is  designed  to  recognize  and  display  solution  modes 
(forced  response,  time  history,  eigensolution )  and  coordinate 
labels  without  user  intervention  beyond  identification  of  co¬ 
ordinate  numbers  selected  for  plotting.  Comparative  and  per¬ 
spective  plots  are  options  (see  Reference  18). 


EXTERNAL  PROGRAN  APPLICATIONS 

Table  7  lists  the  check  cases  used  to  verify  the  capabilities  of 
the  External  Programs.  There  are  six  external  programs  comprising 
the  linked  G400/F389,  Tail  Option  (RIEVA),  modified  E927,  Base 
Program  plot  package,  WABAT  geometry  program/plot  package  (Auto¬ 
matic  Paneling  Technique/Plot  or  APT/Plot),  and  WABAT  potential 
flow  program.  The  check  cases  were  selected  to  verify  key  capa¬ 
bilities,  the  abilities  of  the  programs  to  transmit  data  to  the 
Base  Program  and  to  each  other,  and  to  acquire  data  from  the  Base 
Program,  such  as  plot  data. 
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TABLE  7  -  INPUT  CHECK  CASES  FOR  EXTERNAL  PROGRAMS 
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Coupled  Rotor  Aeroelastic  Analysis  and 
Variable  Inflow  Analysis  G400/F389 


Case  1  verifies  the  ability  of  the  G400  module  to  calculate  an 
articulated  rotor  response.  A  full  scale  rotor  of  22  foot  radius 
was  specified.  Blade  properties  are  typical  of  actual  rotors. 
Uniform  inflow  was  specified.  This  case  uses  specified  input 
control  angles  -  i.e.,  trim  was  not  requested.  Blade  coordinates 
are  flap,  lead,  three  elastic  flatwise,  two  edgewise,  and  one 
torsion  uncoupled  modes. 

Case  2  is  similar  to  Case  1,  except  the  blade  pendulum  vibration 
absorbers  are  specified. 

Case  3  is  an  application  of  G400  which  is  used  to  check  the 
formation  of  the  hub  load  excitation  vector,  rotor  impedance 
matrix,  matrix  1^  containing  the  effects  of  input  HHC  angles  on 
hub  loads  and  blade  stress  matrices  (see  previous  sections  for 
descriptions  of  the  above  terms).  The  data  are  transmitted  to  a 
file.  Base  Program  check  case  13  (Table  6)  uses  the  hub  load 
excitation  vector  and  rotor  impedance  data  to  calculate  system 
response  without  HHC.  Numerical  rotor  trim  is  requested,  as 
opposed  to  analytic  (Bailey  derivative)  trim  (see  Reference  9). 
Analytic  trim  was  invoked  for  the  HHC  ARES  model  applications 
described  subsequently.  In  other  respects,  the  case  is  similar  to 
Case  1 . 

Case  4  is  similar  to  Case  3  except  that  a  teetering  rotor  is 
specified. 

Case  5  exercises  the  F389  program  in  a  stand-alone  mode.  Inflow 
harmonics  data  are  transmitted  to  a  file  which  may  be  read  by  the 
Tail  Option  (RIEVA)  module  to  represent  the  effects  of  the  rotor 
induced  far  wake  at  the  tail. 
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Case  6  exercises  the  F389  program  in  a  stand-alone  mode.  Inflow 
harmonics  data  are  transmitted  to  a  file  by  F389  which  may  be  read 
by  G400 ,  employed  in  a  stand-alone  mode. 

Case  7  is  a  coupled  G400/F389  run.  Two  cycles  of  calculation 
between  G400  and  F389  are  performed  -  i.e.,  the  execution  sequence 
is  G400/F389/G400 ,  performed  twice  (see  Figure  25).  The  rotor  hub 
load  excitation  vector,  rotor  impedance  matrix,  and  HHC  matrix  and 
stress  matrices  are  not  formed.  The  check  of  the  formation  of 
these  matrices  with  variable  inflow  obtained  from  the  coupled 
G400/F389  program  is  obtained  in  the  ARES  model  application 
described  subsequently. 

Tail  Option  Program  (RIEVA) 

Case  1  verifies  the  calculation  of  harmonics  of  tail  loads  when 
rotor  induced  far  wake  data  derive  from  a  file  to  which  inflow 
harmonics  are  transmitted  by  F389.  Figure  32  shows  the  empennage 
airload  variation  with  time. 

Case  2  is  similar  to  Case  1  except  that  the  effects  of  the  far 
wake  are  not  included. 

Rotor  Aeroelastic  Analysis  E927 

Case  1  has  the  rotor  attributes  described  under  Base  Program  check 
Case  12.  Rotor  M,  C,  and  K  matrices  were  transmitted  to  the  Base 
Program . 

Base  Program  Plot  Package 

Case  1  contains  file  data  transmitted  to  a  plot  file  by  the  Base 
Program.  The  file  data  were  calculated  for  the  six  degree-of- 
freedom  isolator  system  described  as  Base  Program  check  Case  11. 
Figures  15,  16,  19,  and  20  illustrate  plots  that  may  be  obtained 
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Figure  32.  Stabilizer  Airloads  Time  History 
in  the  RIEVA  Analysis. 


(see  Reference  18  for  plot  variants). 


WABAT  Geometry/Plot  Program  (APT/Plot) 

Case  1  checks  the  formation  of  panel  geometry  in  the  Automated 
Paneling  Technique  (APT/Plot)  program  from  a  hypothetical  fuse¬ 
lage.  The  fuselage  consists  of  three  sections,  each  made  up  of 
several  panels.  The  sections  consist  of  a  nose  cone,  center 
section  cylinder,  and  reversed  nose  cone,  which  completes  the 
fuselage.  The  APT/Plot  program  assembles  the  three  sections  into 
a  body  resembling  a  fuselage.  The  resulting  geometry  data  are 
stored  for  input  to  the  WABAT  potential  flow  program  or  the  plot 
module  of  APT/Plot. 

Case  2  employs  stored  geometry  data  to  check  the  plot  module  of 
the  APT/plot  program.  The  plot  file  data  contain  panel  geometry 
data  for  a  typical  fuselage.  Plotted  results  are  shown  in  Figures 
27  to  29.  The  results  show  the  capability  of  the  program  to 
permit  viewing  of  the  geometry  from  any  perspective,  and  the 
capability  of  the  program  to  magnify  regions  of  interest. 

Case  3  employs  WABAT  potential  flow  calculations  to  display  flow 
directions  (arrows  in  Figures  27  to  29)  on  the  surfaces  of  fuse¬ 
lage  panels.  It  was  necessary  to  run  the  WABAT  potential  flow 
module  prior  to  obtaining  the  fuselage  flow  plots.  The  case 
illustrates  the  surface  flow  plotting  function  of  the  APT/Plot 
package . 

WABAT  Potential  Flow  Program 

Case  1  uses  APT/Plot  file  data  to  obtain  the  potential  flow 
solution  for  the  fuselage  shown  in  Figures  27  to  29.  Harmonics  of 
fuselage  inflow  at  the  rotor  are  calculated  and  can  be  transmitted 
to  F389 . 
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USER  DOCUMENTATION 


The  following  section  lists  the  documents  which  must  be  used  to 
apply  the  programs  contained  in  the  Coupled  Rotor/Airframe  Vibra¬ 
tion  Analysis  package.  A  brief  description  is  given  of  the 
purpose  and  contents  of  the  documents . 

1.  Cassarino,  S . and  Sopher,  R.:  Coupled  Rotor/Airframe  Vi¬ 
bration  Analysis  Program  Manual.  Volume  I.  User's  and 
Programmer's  Instructions  NASA  CR-165891,  1982. 

The  manual  contains  user  instructions  for  the  Base  Program  des¬ 
cribing  input  preparation  and  output  interpretation.  Input/output 
files  are  described  for  transmitting  data  to  External  Programs. 
Program  structure  and  subroutine  functions  are  described.  Appen¬ 
dixes  describe  the  structure,  file  functions,  and  usage  of  the 
coupled  G400/F389  program,  and  the  E927  program.  Descriptions  are 
given  of  procedures  for  modifying  the  Base  Program.  Sample  test 
cases  are  listed  and  described. 

2.  Volume  2  of  the  above  contains  line  printer  input  and  output 
samples  for  the  Base  Program,  and  coupled  G400/F389  and  E927 
programs . 

3.  Bielawa,  Richard,  L.  :  Aeroelastic  Analysis  for  Helicopter 
Rotors  with  Blade  Appended  Pendulum  Vibration  Absorbers  - 
Mathematical  Derivations  and  Program  User's  Manual.  NASA 
CR-165896, 1982 . 

The  mathematical  basis  of  modifications  to  the  G400  analysis 
comprising  teetering  rotor  and  pendulum  absorber  modifications  and 
the  method  used  to  calculate  the  rotor  impedance  are  described. 
The  complete  mathematical  basis  of  the  analysis  prior  to  modifica¬ 
tions  is  described  in  Reference  8  listed  also  as  item  4  below. 
Program  hierarchy  and  subroutine  functions  are  described.  Input 


preparation  is  described. 


4.  Bielawa,  R.  L.  :  Aeroelastic  Analysis  for  Helicopter  Rotor 
Blades  with  Time-Variable,  Nonlinear  Structural  Twist  and 
Multiple  Structural  Redundancy  -  Mathematical  Derivation  and 
Program  User's  Manual.  NASA  CR-2638,  1976. 

5.  Egolf,  T.  A.,  and  Landgrebe,  A.  J.:  A  Prescribed  Wake  Rotor 
Inflow  and  Flow  Field  Prediction  Analysis.  NASA  CR-165894, 
1982. 

This  manual  describes  the  mathematical  basis  of  the  rotor  induced 
variable  inflow  analysis,  F389,  functions  of  program  subroutines, 
and  input  preparation.  The  analysis  is  a  single  rotor  version  of 
the  analysis  listed  below  as  item  6,  which  provides  further 
details  of  the  mathematical  basis  of  the  method. 

6.  Landgrebe,  A.  J.,  and  Egolf,  T.  A.:  Rotorcraft  Wake  Analy¬ 
sis  for  Prediction  of  Induced  Velocities.  USAAMRDL  TR 
75-45,  1976. 

7.  Cassarino,  Sebastian,  J.  ,  and  Mouzakis,  Timoleon:  Bifilar 

Analysis  User's  Manual  -  Volume  II.  NASA  CR-159228,  1980. 

The  manual  describes  input  preparation,  program  structure,  and 
functions  of  subroutines  of  the  E927  program  including  a  bifilar 
vibration  absorber.  Modifications  were  made  to  the  bifilar 
version  of  E927  program  for  the  present  contract  to  eliminate  the 
bifilar  capability  and  to  configure  the  E927  rotor  as  a  sub¬ 
structure  for  use  by  the  Base  Program.  Appendix  B  of  item  1  above 
describes  changes  users  must  make  to  E927  input  to  employ  the 
modified  program  and  specifies  the  file  for  transmitting  data  to 
the  Base  Program.  Sample  inputs  and  outputs  are  provided  in  item 
1  above.  The  mathematical  basis  of  rotor  aeroelastic  analysis 
E927  is  described  in  Reference  6,  listed  as  item  8  below. 
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8.  Johnston,  R.  A.,  and  Cassarino,  S.  J.;  Helicopter  Rotor 

Stability  Analysis,  USAAMRDL  TR-75-40,  U.S.  Army,  January 
1976  . 

9.  Gangwani ,  Santu  T. :  A  Doublet  Lattice  Method  for  the 

Determination  of  Rotor  Induced  Empennage  Vibration  Airloads 

Analysis  Description  and  Program  Documentation.  NASA 
CR-165893 ,  1982. 

The  manual  describes  the  basis,  program  structure,  subroutine 
functions,  and  input  preparation  for  use  of  the  Rotor  Induced 

Empennage  Vibration  Airloads  (RIEVA)  program.  Applications  are 
described  and  discussed.  The  file  for  accepting  variable  inflow 
data  from  the  F389  program  is  specified  in  item  1  above. 

10.  Studwell,  R.  E.:  User's  Manual  for  the  Coupled  Rotor/Fuse¬ 

lage  Vibration  Analysis  Graphics  Package.  NASA  CR-165897, 
1982  . 

The  manual  describes  the  user  responses  to  plot  package  messages 
for  plotting  results  transmitted  by  the  Base  Program.  Plot  option 
results  are  illustrated.  Program  function  and  subroutine  struc¬ 
ture  are  described.  The  plot  data  transmission  file  format  is 
specified  in  the  text  of  item  1  above. 

11.  Studwell,  R.  E. :  User's  Manual  for  the  Automated  Paneling 

Technique  (APT)  and  the  Wing  Body  Aerodynamic  Technique 
(WABAT)  Programs.  NASA  CR-165895,  1982. 

This  manual  describes  the  preparation  of  input  for  the  WABAT 
geometry  package  (APT/Plot)  and  the  WABAT  potential  flow  program. 
The  manual  describes  the  user  responses  to  plot  fuselage  geometry 
and  fuselage  surface  flows.  Sample  problems  are  illustrated. 
Input  data  requirements  for  the  geometry  package  (APT/Plot)  are 
described.  The  basis  of  the  automated  paneling  technique  is 
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described  in  Appendix  A  of  Reference  13  (see  item  12  below).  The 
potential  flow  analysis  basis  is  summarized  in  Appendix  B  of 
Reference  13  (see  item  12  below). 

12.  Sheehy,  T.  W. :  A  Method  for  Predicting  Helicopter  Hub  Drag. 

( USAAMRDL  Contract  DAAJ02-74-C-0050 . )  USAAMRDL-TR-75-48 , 
1975  . 


CORRELATION  OF  THEORY  AND  TEST 


In  addition  to  development  of  the  Coupled  Rotor/Air  frame  Vibration 
Analysis,  this  contractual  effort  included  a  limited  correlation 
study.  The  correlation  study  involved  using  wind  tunnel  data  from 
a  dynamically  scaled  UH-60A  rotor  system  coupled  with  the  NAL.. 
Langley  aeroelastic  rotor  experimental  system  (ARES)  powered  test 
stand.  Tests  were  conducted  with  and  without  higher  harmonic 
control  ( HHC )  employed.  Model  balance  loads,  rotor  bending 
moments  and  vibrations  in  the  ARES  structure  were  measured.  A 
shake  test  was  performed  to  obtain  the  modal  attributes  of  the 
system  for  input  to  the  Coupled  Rotor/Airframe  Vibration  Analysis. 

Shake  Test  of  the  ARES  Model 

A  shake  test  was  conducted  in  the  NASA  Langley  Transonic  Dynamics 
Tunnel  (TDT)  to  determine  the  support  system  modal  properties 
required  for  vibration  analysis  predictions.  The  support  system 
was  provided  by  the  Langley  Research  Center  Aeroelastic  Rotor 
Experimental  System  (ARES)  (see  Reference  19).  The  test  procedure 
consisted  of  exciting  the  model  by  an  electromagnetic  shaker  in 
several  different  directions.  Random  excitations  were  applied  to 
the  model.  Data  were  recorded  on  14  track  FM  magnetic  tapes  and 
consisted  of  outputs  from  eight  hub  accelerometers,  three  body 
accelerometers,  six  rotor  balance  channels,  and  the  force  load 
cell . 
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The  shake  test  data  were  processed  using  a  dual  channel  digital 
analyzer  (HP  5420)  to  obtain  transfer  functions  -  e.g.  accel¬ 
eration/applied  force  harmonics  -  of  the  measurements  with  respect 
to  the  applied  forces.  The  ARES  modal  attributes  presented  in 
this  section  and  those  used  in  the  correlation  study  are  based  on 
the  random  excitations.  Table  8  lists  the  characteristics  of  nine 
ARES  modes  covering  a  frequency  range  from  21  to  73  Hz.  The 
frequencies  of  interest  are  located  near  45  Hz  which  is  4/rev.  for 
the  four  bladed  model  rotor.  The  natural  frequency,  generalized 
macs,  percent  critical  damping,  six-component  mode  shape  at  the 
hub,  and  three  translational  components  of  modal  response  at  the 
airframe  accelerometer  location  for  each  mode  are  listed  in  Table 
8. 

Modes  1,  2,  3,  7,  8  and  9  were  obtained  from  a  vertical  off  center 
excitation  of  the  hub.  A  centered  lateral  excitation  run  defined 
the  properties  for  modes  4  and  5.  A  centered  longitudinal  excita¬ 
tion  run  defined  mode  6,  The  generalized  masses  for  each  of  these 
nine  modes  were  determined  from  single  degree  of  freedom  analyti¬ 
cal  models  of  the  system.  The  damping  in  each  mode  was  obtained 
by  peak  fit  using  the  HP  5420.  Imaginary  components  of  the 
transfer  functions  were  used  to  define  the  mode  shapes  at  the 
rotor  hub  and  airframe. 

Wind  Tunnel  Test  Program 

The  test  program  was  conducted  in  the  NASA  Langley  Transonic 
Dynamic  Tunnel  (TDT).  The  model  rotor  was  install  1  on  the 
Langley  Research  Center  Aeroelastic  Rotor  Experimental  System 
(ARES),  Reference  19.  Figure  33  shows  the  ARES  model  with  the 
test  rotor. 

A  four-bladed  articulated  rotor  (1/6  scale)  was  installed  on  the 
ARES  model.  The  dynamically  scaled  blades  were  based  on  UH-60A 
blades.  The  same  blades  were  used  for  baseline  studies  in 
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Reference  20  which  contains  details  of  the  blade  properties.  The 
blades  were  scaled  for  operation  at  full-scale  Mach  numbers  in 
Freon.  Model  Reynolds  numbers  were  40  percent  of  full-scale 
values.  The  blades  have  a  56.22-inch  radius,  a  3.625  inch  chord, 
and  a  -12  degree  twist.  The  blades  have  an  SC1095  airfoil  section 
i’iboard  of  the  50  percent  and  outboard  of  the  85-percent  radius 
stations.  The  SC1095R8  airfoil  is  used  between  the  50-  and 
35-percent  span  stations. 

To  achieve  vibration  reductions  with  HHC,  various  types  of  higher 
harmonic  controllers  were  tested.  The  control  algorithms  are 
described  in  References  11  and  21.  Comparisons  of  theoretical 
predictions  obtained  with  the  Coupled  Rotor/Airframe  Vibration 
Analysis  were  made  with  test  data  for  the  adaptive  closed-loop 
controller  listed  as  controller  2  in  Table  1  of  Reference  21.  In 
this  controller,  the  T  matrix  was  identified  and  updated  on-line 
and  "cautious"  terms  were  provided  (see  Reference  21). 

Application  of  the  Analysis 

The  Base  Program  was  applied  in  its  forced  response  solution  mode, 
FR1 .  The  forcing  frequency  used  in  this  solution  mode  for  the 
present  application  was  44.7  Hz  which  is  the  4/rev.  frequency  for 
the  four  bladed  rotor  model  at  its  operating  speed  of  670  RPM. 

Rotor  Representation.  The  rotor  effects  were  represented  by 
substructure  RE3  which  models  a  multibladed  aeroelastic  rotor  in 
forward  flight  or  hover.  This  component  reads  the  rotor  hub 
excitation  force  and  moment  vector  and  the  rotor  impedance  matrix 
and  HHC  and  blade  stress  matrices  computed  by  rotor  aeroelastic 
analysis  G400.  In  this  application  G400  was  executed  iteratively 
with  External  Program  F389  to  account  for  rotor  induced  variable 
inflow.  The  iterative  sequence  was  G400-F389-G400  -  i.e.,  one 
coupled  cycle  of  these  programs  was  run  (see  Figure  25). 


127 


The  rotor  was  analytically  trimmed  to  achieve  the  test  values  of 
lift  and  shaft  angles  for  advance  ratios  from  0.20  to  0.40.  Blade 
properties  were  taken  from  Reference  20.  The  resulting  blade 
frequencies  are  listed  in  Table  9.  Two  flatwise,  one  edgewise, 
and  one  torsion  mode(s)  were  employed  in  the  computations. 
Full-scale  Reynolds  number  data  for  the  SC1095  and  SC  1095R8 
airfoils  were  used  in  rotor  aeroelastic  analysis  G400.  The  thrust 
coefficient  CT/o  was  equal  to  0.07  for  all  cases. 

Rotor  Impedance  Sensitivity  to  Perturbation  Amplitudes.  -  Ex¬ 
ternal  Program  G400  obtains  the  rotor  impedances  with  respect  to 
hub  motions  and  HHC  angles.  The  perturbation  amplitudes  were 
selected  to  achieve  accurate  values  of  the  numerically  obtained 
impedance  elements.  Applications  showed  that  there  was  a  practi¬ 
cal  limit  on  the  size  of  the  perturbation  amplitudes  for  accurate 
results.  Accurate  results  were  verified  by  comparisons  of  G400 
numerical  results  with  analytically  derived  results  for  rotor 
impedance  and  HHC  control  matrix  elements. 

Figures  34  and  35  show  sensitivity  of  the  z-...(l,l)  and  h  (3,3) 

-L  .L  1  in 

elements  of  the  rotor  impedance  and  h^  matrices  to  perturbation 
amplitude.  Based  on  these  figures,  the  perturbation  amplitudes 
selected  were  152  m/sec2  (500  ft/sec2)  for  hub  linear  acceler¬ 
ation,  50  rad/sec2  for  hub  angular  acceleration,  and  0.10  degree 
for  HHC  angles. 

Support  System  Representation.  -  The  ARES  modal  attributes  listed 
in  Table  8  were  incorporated  into  the  analysis  by  the  use  of  modal 
substructure  MSI.  Each  ARES  mode  was  represented  by  one  MSI 
substructure.  The  connection  point  between  the  rotor  hub  and 
support  system  was  defined  as  the  first  node.  The  airframe 
location  where  vibrations  were  measured  was  taken  as  the  second 
node.  The  airframe  node  required  the  use  of  a  generalized  force 
GF1  type  of  substructure. 
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TABLE  9  -  CALCULATED  BLADE  NATURAL  FREQUENCIES  FOR 
THE  MODEL  ROTOR  (PER  REV) 


Mode  No . 

1 

2 

3 

Flatwise 

2 . 7 

4 . 7 

7 . 7 

Edgewise 

5 . 0 

14.2 

Torsion 

6 . 0 

14.9 

TABLE  10  -  COMPARISON  OF  ARES  AIRFRAME  AND  MODEL  ROTOR 
BLADE  NATURAL  FREQUENCIES 


ARES  SHAKE  TEST  DATA 
(REFERRED  TO  FIXED  SYSTEM) 


CALCULATED 
BLADE  MODES 

(REFERRED  TO  MOVING  SYSTEM) 


No . 

Description 

io /Q 

Ei 

- 

1 . 9 

- 

2.3 

- 

2 . 5 

4 

Lateral  Roll 

3 . 4 

5 

Lateral  Roll 

4.4 

6 

Longitudinal  Pitch 

5.2 

Description 

Flatwise  1 

2.7 

Flatwise  2 

4.7 

Edgewise  1 

5.0 

3  Mode 
Solution 
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LONGITUDINAL  PERTUBATION 


0  500  1000 

| - 1 - 1 - \ 

0  100  200  300 

Figure  34.  Sensitivity  of  Rotor  Impedance 
to  Hub  Perturbation  Amplitude. 
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Specified  Higher  Harmonic  Control.  -  Four  values  of  pitch  ampli¬ 
tude  were  specified  to  determine  the  effects  of  non-optimal 
control  angle  inputs  on  vibrations  and  loads.  These  four  cases 
were  obtained  by  specifying  e3c  -  +  1°,  and  -1°,  and  e3g  =  +  1°, 
and  -1°.  These  specified  inputs  are  non-optimal  because  they  are 
not  determined  from  the  requirement  that  they  must  minimize 
vibrations . 

Optimal  Higher  Harmonic  Control.  -  The  analytical  optimal  control 

state  is  derived  from  equation  (93),  which  is  essentially  the  same 

as  the  algorithm  presented  in  reference  11.  The  control  angle 

vector  has  six  components  listed  in  equations  (78)  and  79)  -  i.e. 

cosine  and  sine  components  of  3-,  4-,  and  5-per-rev.  inputs. 

Three  airframe  vibrations  (longitudinal,  lateral,  and  vertical) 

were  selected  for  inclusion  in  the  performance  index,  J,  to  be 

minimized  (see  equation  (90)).  Elements  of  the  weighting  matrix 

w  were  20,000,  20,000,  and  150,000  -  i.e.,  vertical  vibration 

z 

reduction  was  given  more  weight  than  rotor  inplane  vibration 
reductions.  Diagonal  elements  of  the  control  weight  matrix  used 
to  limit  the  amplitudes  of  HHC  control  inputs  were  each  equal  to 
.01.  Analytical  control  inputs  were  found  to  be  less  than  1.5° 
for  all  conditions  except  m  =0.5  where  03s  was  predicted  equal  to 
3°.  The  test  HHC  inputs  were  limited  to  a  maximum  of  1.5°.  The 
analytical  controller  employed  an  analytical  determination  of  the 
T  matrix  (see  Higher  Harmonic  Control  for  Vibration  Reduction). 
Noise  was  not  represented. 

Correlations  for  Baseline  Conditions 

Airframe  Vibrations.-  Comparison  of  predicted  vibration  levels  for 
a  three  mode  and  nine  mode  solution  versus  ARES  test  data  are 
shown  in  Figure  36.  Both  solutions  adequately  predict  the  trend 
of  increasing  vibration  levels  with  increasing  advance  ratio.  The 
modes  used  for  the  three  mode  solution  are  those  modes  with 
frequencies  closest  to  the  rotor  4/rev.  which  should  be  the  primary 
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contributions  to  vibration  levels.  The  nine  mode  solution  was 
employed  in  an  attempt  to  improve  the  correlation  with  test  data 
and  assess  the  sensitivity  to  modal  inputs.  As  would  be  expected, 
the  additional  modes  resulted  in  higher  absolute  values  of  vibra¬ 
tion  levels,  however,  despite  significant  fuselage  modal  partici¬ 
pation  of  modes  2,  3,  8,  and  9,  the  large  contribution  to  vibra¬ 
tion  levels  was  not  expected.  In  fact,  inclusion  of  the  addi¬ 
tional  six  modes  results  in  predicted  vibration  levels  signifi¬ 
cantly  higher  than  measured  values.  All  modes  used  were  obtained 
from  random  shake  test  excitations  which  historically  have  proved 
less  accurate  than  swept  sine  results.  Since  the  vibration 
predictions  appear  to  be  very  sensitive  to  modal  characteristics, 
a  higher  level  of  accuracy  in  representing  the  fixed  system  modes 
is  desirable. 

Rotor  Blade  Moments.  -  Blade  vibratory  bending  and  torsion  moment 
trends  with  airspeed  predicted  by  the  analysis  (see  Figure  37) 
were  examined.  Edgewise  moments  are  overpredicted  whereas  the 
torsion  moments  are  underpredicted  by  more  than  50  percent.  Blade 
moment  variations  with  radius  shown  in  figure  38  for  (j=  0.30  also 
have  the  same  behavior . 

The  torsion  test  data  in  Figure  38  appear  to  contain  errors.  Test 
data  for  C^/o  =  .075  and  (j  =  0.3  derived  from  Reference  20  are 
plotted  in  Figure  39.  The  test  flatwise  and  edgewise  moments  show 
a  slight  increase  in  comparison  with  the  curves  in  Figure  38.  The 
torsion  correlation  is  good,  suggesting  that  the  torsion  test  data 
in  Figure  38  contain  errors. 

Figures  40  and  41  show  the  harmonics  of  the  blade  moments  at  mid 
span  and  )j=  0.30.  Figure  41  contains  test  results  from  Refer¬ 
ence  20.  The  flatwise  moment  harmonics  are  predicted  adequately 
by  the  analysis. 
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Figure  38.  Radial  Distributions  of  Blade 
Moments  at  0.3  Advance  Ratio. 


Ratio  with  Test  Data  from  Reference  20. 


The  edgewise  harmonics  correlate  satisfactorily  except  at  5-  and 
6-per-rev.  The  analysis  predicts  a  large  5/rev.  moment  and  a 
small  6/rev.  moment.  The  large  5/rev.  component  accounts  for  the 
overprediction  of  the  edgewise  moments,  (see  Figures  37  and  39). 
This  5/rev.  amplification  is  probably  due  to  resonance  of  this 
mode  with  4/rev.  rotor  hub  excitation,  as  explained  below. 

Pitch  link  load  comparisons  were  similar  to  the  torsional  moment 
curves  in  Figure  37. 

The  conclusion  is  that  the  blade  moment  trends  are  adequately 
predicted  by  the  analysis.  Absolute  values  appear  to  depend  on 
how  accurately  the  blade  parameters  are  modelled,  especially, 
blade  frequencies. 

Hub/Balance  Force  Correlations.  -  Analytical  4-per-rev.  hub  force 
and  balance  test  force  trends  with  airspeed  are  shown  in  Figure 
42.  A  direct  comparison  of  hub  loads  is  not  made  because  the  test 
data  consist  of  balance  loads  while  the  analysis  data  consist  of 
hub  loads.  A  dynamic  calibration  matrix  is  required  to  relate  the 
4-per-rev.  loads  at  the  model  hub  to  the  4-per-rev.  loads  sensed 
by  the  balance.  The  calibration  matrix  was  not  available  in  this 
correlation  study.  Figure  42  also  shows  that  the  longitudinal  hub 
force  from  analysis  has  a  greater  sensitivity  to  airspeed  than  the 
longitudinal  balance  force.  The  lateral  and  vertical  hub  force 
trends  with  airspeed  correlate  fairly  well  with  the  corresponding 
balance  loads.  Since  no  direct  comparisons  of  hub  loads  are 
available,  no  conclusions  can  be  made  on  the  accuracy  of  the 
analysis,  except  that  the  analysis  appears  to  show  the  proper 
trend  of  hub  loads  with  airspeed. 

Effect  of  Hub  Displacements.  -  An  area  of  interest  is  the  sensi¬ 
tivity  of  blade  moments  and  hub  loads  to  support  system  flexi¬ 
bility.  If  the  support  system  modal  attributes  are  such  that 
blade  loads  are  substantially  different  from  those  with  a  fixed 
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Figure  42.  Comparison  of  Variations  with  Advance  Ratio  of 
Predicted  Hub  Loads  and  Measured  Balance  Loads 


hub,  then  the  blade  analysis  should  properly  account  for  effects 
of  hub  displacements. 

Figure  43  shows  that  blade  edgewise  moments  are  affected  by  hub 
flexibility,  and  that  hub  motion  increases  the  alternating  edge¬ 
wise  moment  at  mid  span  by  50  percent.  Table  10  reveals  the 
source  of  the  increase  in  edgewise  moments.  The  blade  has  an 
edgewise  mode  at  5/rev.  which  transforms  into  4/rev.  and  6/rev.  in 
the  fixed  system.  Only  the  4/rev.  component  is  transmitted  to  the 
fixed  system  for  this  4-bladed  rotor  system.  This  transformed 
4/rev.  edgewise  mode  is  at  resonance  with  the  4/rev.  rotor  excita¬ 
tion  and  may  couple  with  the  fixed  system  lateral  roll  mode  at 
4.4/rev.  This  speculation  is  reinforced  by  Figure  44  which  shows 
that  the  in-plane  hub  forces  are  more  sensitive  to  hub  displace¬ 
ment  than  the  vertical  hub  force . 

The  deduction  is  that  hub  participation  can  have  a  substantial 
effect  on  blade  moments  and  hub  loads.  The  accurate  prediction  of 
hub  loads  is  sensitive  to  accurate  characterization  of  support 
system  modal  attributes  and  rotor  blade  frequencies.  The  lack  of 
agreement  of  edgewise  moments  with  test  data  is  likely  to  be  due 
to  inaccurate  characterization  or  rotor  blade  natural  frequencies 
in  the  analysis.  Because  rotor  loads  are  the  primary  excitation 
sources  of  coupled  rotor/airframe  systems,  these  sensitivities 
show  the  importance  of  accurate  modal  characteristization  of  both 
the  rotor  and  airframe  for  vibration  prediction. 

Correlation  with  HHC  Non-optimal  Inputs 

The  Coupled  Rotor/Airframe  Vibration  analysis  was  executed  to 
predict  airframe  vibrations,  rotor  bending  moments  and  hub  loads 
utilizing  arbitrary  higher  harmonic  control  inputs.  Comparison 
with  wind  tunnel  measurements  for  the  same  inputs  is  discussed 
below. 
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Airframe  Vibrations.  -  Figure  45  shows  the  effect  of  3/rev.  HHC 
pitch  inputs  on  airframe  4/rev.  vibration.  Also  shown  is  the 
analytical  estimate  of  vibration  obtained  from  the  Coupled  Rotor/ 
Airframe  Analysis  with  HHC  and  a  three  mode  representation. 
Absolute  values  are  underpredicted  and  the  trend  in  vibration 
reduction  with  HHC  application  is  not  predicted  uniformly.  The 
nine  mode  solution  results  shown  in  Figure  46  in  general  agree 
better  with  the  test  data  and  shows  uniformly  higher  vibration 
levels  than  the  three  mode  solution. 

Figure  47  shows  vibration  vector  corrrelations  for  the  airframe 
4/rev.  longitudinal  vibrations.  For  all  four  HHC  input  values, 
the  analysis  is  able  to  predict  the  phase  changes  due  to  HHC 
inputs.  Amplitude  comparisons  are  the  same  as  in  Figure  45. 

Rotor  Blade  Moments.  -  Figure  48  shows  the  effect  of  HHC  inputs  on 
blade  bending  moments.  This  figure  indicates  that  bending 
moments,  especially  edgewise  moments,  are  sensitive  to  the  HHC 
inputs.  Harmonic  analysis  of  the  moments  (see  Figure  49)  indi¬ 
cates  the  differences  are  primarily  due  to  large  3/rev.  flatwise 
and  6/rev.  edgewise  contributions  for  the  HHC  applications.  The 
predicted  blade  moment  distribution  for  baseline  and  HHC  input 
case  are  provided  in  Figure  50.  The  sensitivity  to  HHC  inputs 
observed  in  the  measured  data  is  not  predicted. 

The  sensitivity  of  the  measured  blade  loads  to  a  3/rev.  input  may 
be  due  to  impurities  in  the  3/rev.  pitch  input  and  the  proximity 
of  the  test  blade  edgewise  mode  to  6/rev.  The  impurities  are 
assumed  to  give  rise  to  a  6/rev.  excitation  which  excites  the 
6/rev.  edgewise  mode.  This  observation  was  made  in  Reference  11. 
The  analysis  which  models  a  3/rev.  HHC  input  as  a  3/rev.  pitch 
change  is  not  able  to  predict  this  effect. 
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Effects  of  Specified  HHC  on  Amplitudes  of 
Airframe  Vibrations  -  Three  Mode  Solution. 
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Figure  47.  Effects  of  Specified  HHC  on  Longitudinal 
Vibration  Vectors  -  Three  Mode  Solutions, 
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Figure  48.  Effects  of  Specified  HHC  on  Measured 
Blade  Moment  Distributions. 
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Figure  49.  Effects  of  Specified  HHC  on  Harmonics 
of  Measured  Blade  Moments. 
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Figure  50.  Effects  of  Specified  HHC  on  Predicted 
Blade  Moment  Distributions. 
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Hub/Balance  Force  Correlations.  -  Figure  51  shows  side-by-side 
comparisons  of  predicted  hub  forces  with  test  balance  forces . 
Variations  with  HHC  inputs  are  not  well  correlated.  Because  no 
direct  hub  load  comparisons  are  made,  no  conclusions  are  drawn. 

Optimal  HKC  Solutions 

Analytical  and  test  results  are  compared  for  the  effects  of 
optimal  values  of  HHC  inputs  on  airframe  vibrations,  blade 
moments,  and  hub  loads. 

Airframe  Vibration.  -  Predicted  and  measured  airframe  vibrations 
for  the  baseline  model  and  optimal  HHC  inputs  as  functions  of 
advance  ratio  are  shown  in  Figure  52 .  The  analysis  predicts  the 
proper  trends  with  HHC  application  and  advance  ratio  increases  for 
both  the  vertical  and  longitudinal  accelerations.  The  test  data 
indicate  that  the  HHC  inputs  are  not  effective  in  reducing  lateral 
vibrations  in  contrast  to  the  analytical  predictions. 

Blade  Bending  Moments.  -  Both  analysis  and  test  show  that  for  the 
conditions  selected,  the  spanwise  blade  bending  moment  distribu¬ 
tions  are  not  sensitive  to  HHC  inputs.  Figure  53.  The  disagree¬ 
ment  in  the  torsion  moments  is  probably  due  to  errors  in  the  test 
data  as  noted  in  the  correlation  of  baseline  conditions. 

Figure  54  shows  that  predicted  flatwise  and  edgewise  moments  are 
consistent  with  measured  values  and  that  the  trends  with  advance 
ratio  are  predicted  correctly.  The  trend  for  torsional  moments  is 
consistent  but  the  predicted  magnitudes  are  consistently  lower 
than  measured. 

Figure  55  shows  comparisons  between  test  results  and  predicted 
blade  moment  harmonics  at  midspan  and  p=  0.30.  The  torsion 
moments  will  not  be  discussed  because  the  test  data  probably 
contain  errors.  The  flatwise  and  edgewise  harmonics  of  moments 
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Figure  52.  Effects  of  Optimal  HHC  on  Fuselage  Vibrations. 
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Figure  56.  Effects  of  Optimal  HHC  on  Comparisons  of 
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correlate  fairly  well.  The  overprediction  of  the  5/rev.  and 
underprediction  of  the  6/rev.  edgewise  components  may  be  due  to 
the  input  blade  edgewise  mode  frequency  being  5/rev.  while  the 
actual  frequency  is  closer  to  6/rev. ,  as  discussed  in  the  pre¬ 
ceding  text. 

Hub/Balance  Force  Comparisons.  -  Figure  56  shows  that  the  analysis 
predicts  a  large  effect  of  the  HHC  inputs  on  the  longitudinal  hub 
force  which  is  not  evidenced  in  the  measured  balance  loads.  The 
effects  on  lateral  and  vertical  forces  are  more  consistent  but  the 
measured  balance  vertical  forces  are  more  sensitive  to  HHC  inputs 
than  the  predicted  effect  on  vertical  hub  forces. 

CONCLUDING  REMARKS 


1.  A  Coupled  Rotor/Airframe  Vibration  Analysis  has  been 
developed  to  serve  as  a  preliminary  design  tool  for  pre¬ 
dicting  helicopter  vibrations  and  a  research  tool  to 
quantify  the  effects  of  structural  properties,  aerodynamic 
interactions  and  vibration  reduction  devices  on  vehicle 
vibration  levels. 

2.  The  analysis  developed  employs  a  substructure  method  which 
provides  an  efficient  stand-alone  mode  operation  of  a  Base 
Program  while  retaining  a  high  level  of  aerodynamic  and 
dynamic  sophistication  by  linking  output  of  External 
Programs . 

3.  The  Base  Program  operating  in  a  stand-alone  mode  can 
efficiently  be  used  to  study  design  changes  in  fuselage 
properties  and  vibration  reduction  devices. 
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4. 


The  substructure  basis  of  the  Base  Program  is  well  suited  to 
allowing  the  Base  Program  to  introduce  data  from  user 
provided  External  Programs  of  the  rotor  system.  The  sub¬ 
structure  basis  also  will  facilitate  growth  of  the  Base 
Program  to  include  additional  internal  components. 

5.  The  Rotor/Airframe  coupling  in  the  analysis  is  accomplished 
using  an  impedance  matching  technique. 

6.  The  current  External  Programs  supporting  the  Base  Program 
provide  the  capability  to  model  articulated,  hingeless  and 
teetering  rotor  systems,  higher  harmonic  control  inputs, 
main  rotor  wake/empennage  excitation  sources  and  hover, 
transition  and  forward  flight  regimes. 

7.  A  limited  correlation  study  with  a  dynamically  scaled  wind 
tunnel  model  indicates: 

a.  The  analysis  predicts  the  proper  trends  in  vibration 
levels  as  a  function  of  advance  ratio. 

b.  The  analysis  predicts  the  proper  trends  in  vibration 
levels  and  blade  bending  moments  with  application  of 
higher  harmonic  control. 

c.  Absolute  values  of  predicted  vibration  levels  and  blade 
bending  moments  are  sensitive  to  airframe  modal  char¬ 
acteristics  and  were  not  uniformly  consistent  with 
measured  values. 
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APPENDIXES 


Appendix  A  -  Illustration  of  System  Synthesis 
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Mapping  From  Dependent  Coordinates  To  Independent  Coordinates 


Component  displacements  are  equated  at  connection  node  1 
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Assembled  system  equation  of  motion 
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MXj  +  K  Xj  =  0 

Recovery  of  Component  Responses  and  Reactions 
Components  responses 
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(A33  ) 


Component  reactions 


Rd  =  FD  -  (mdxd  +  W 


( A3  5  ) 
( A36  ) 


Rotation  Transformation 


(Global  Axis) 


(Local  Axis) 


Sketch  ( A2 ) 
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Organization  of  Coordinates  XQ  and  X^ 

Elem  4 
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Composite  Transformation 


For  the  case  where  substructure  and  global  axis  orientations  are 

the  same 

P0  =  I  (A43 ) 

*D  =  P0X0  =  xo 

P  =  P2  ( A45 ) 

=  P2  ( A46 ) 
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Appendix  B  -  Substructure  Equations 


Fixed  Absorber  FA1 

See  Figure  5  for  schematic  of  absorber. 
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Matrices  for  the  fixed  absorber  are  shown  in  Figure  Bl. 

The  force  vector  is  (see  Figures  2  and  5) 

{ f i  =  {O}  (Bl) 


The  reaction  load  vector  (see  Figure  2  for  positive  load  direc¬ 
tions)  is 
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Generalized  Force  GF1 


See  Figure  6  for  schematic  of  substructure. 
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The  matrices  for  the  GF1  substructure  are 


[ml 

= 

[0] 

( B5  ) 

[c] 

= 

[0] 

( B6  ) 

DO 
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( B7  ) 

The  specified  load  vector  at  the  connection  node  is 


{f}  =  {f1,f1,f1,m,,m.,m,} 

*  xl  yl  zl  xl  yl  zlJ 

The  reaction  load  at  the  connection  node  is 


( B8 ) 
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( B9  ) 
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Beam  Model  BM1 


See  Figure  7  for  schematic  of  beam. 
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Figure  B2  illustrates  the  mass  matrix  of  the  system.  Figure  B3 
illustrates  the  stiffness  matrix.  The  applied  force  vector  is 


{f}  =  *0} 
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The  reaction  load  vector 
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Directions  of  force  components  at  a  connection  node  are  shown  in 
Figure  2 . 


Properties  of  the  beam  are 
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Constraint  Substructure  CN1 
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M  x  1 


( B18  ) 


m  =  ( o } 


Constraint  substructure  CN1  is  a  special  type  of  modal  substruc¬ 
ture.  Equation  (20)  is  the  basis  of  the  constraint  procedure.  As 
an  example,  if  it  is  desired  to  constrain  -  i.e.  eliminate  as 
coordinates  -  all  physical  displacements  at  a  connection'  node 
except  two  -  say,  and  6^  -  equation  (20)  takes  the  special  form 


GLOBAL 

ORIENTATION 

AXIS 


(B19) 


This  is  obtained  by  assuming  [E]  =  I  and  replacing  elements  of  the 
modal  matrix  by  zeroes  where  it  is  desired  to  constrain  coord¬ 
inates,  and  ones  where  it  is  desired  to  retain  coordinates.  The 
modal  coordinates  q^  and  q2  are  identical  to  v^,  02.  Analogous 

transformations  apply  to  other  selection  combinations. 

Aeroelastic  Rotor  Type  2  ( RE2 ) 

Rotor  coordinates  are  listed  in  Table  4.  Equations  of  motion  are 
defined  in  reference  6. 
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Modal  Substructure  MSI 
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X 
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Matrices 
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( B20  ) 

[c]  - 

[°n] 
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[k]  = 

Ikul 

( B22  ) 

{f}  = 

(0) 

( B23  ) 

At  a  connection  node,  the  generalized  reaction  load  vector  is 

T 

{ Q)r  =  (il  {r}  ( B24 ) 


The  modal  vector  at  a  connection  node  is 


{ii  *  <Yxi)r  ( y  yi ^ l '  *yyl^l'  ^yei^l'  ^ y02 ^1 ^ ^ y03 ^1 ! 


( B25  ) 


The  elements  of  the  modal  vector  are  displacements  which  have  the 
positive  directions  shown  in  Figure  2. 


The  reaction  load  vector  at  a  connection  node  is 


{r} 


{rxl'  ryl '  rzl'  rmxl '  rmyl '  rmzl^ 

The  displacement  vector  at  a  connection  node  is 


{x}  =  {*}  {q} 

T 

{x}  =  {  u^,  v^,  w^,  ei,  e2 ,  03)  local  orientation 

1  x  "  axes 


( B26 ) 


( B27 ) 
( B28  ) 


If  a  physical  substructure  is  not  connected  to  a  node  of  the  modal 
substructure,  obtain  the  physical  displacements  of  the  node  by 
connecting  a  GF1  substructure  to  the  node  with  zero  applied  loads 
specified  in  GF1 . 
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Horizontal  Linear  Bifilar  BF1 
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The  multi-blade  coordinates  y  and  y  are 

c  s 

bifilar  arm  lead  physical  displacements  by 


Y  =  (Yc  c°s<ii  +  y  s  sinijj )  ( 2/N ) 

See  Figures  8  and  2  for  positive  directions  of  y 
node  displacements.  Figure  B4  shows  the  non-zero 
the  mass,  damping  and  stiffness  matrices. 


( B29 ) 


( B30 ) 


related  to  the 

(B31) 

and  connection 
submatrices  of 
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The  applied  load  vector  is 

{f}  =  {0}  (B32) 

The  reaction  load  vector  is 


{ r }  ={Of  O,  r,,r,,r,,r  ,,r  ,,r  ,  } 

11  1  xl  yl  zl  mxl  myl  mzl' 


( B33  ) 


Figure  2  illustrates  the  positive  directions  of  the  reaction 


loads . 


Properties 

m,  -  bifilar  arm  mass 
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R  -  equivalent  hinge  offset 

r  -  equivalent  pendulum  length 

n  -  tuning  ratio 

=  (R/r)H 
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Q  -  rotation  speed 

N  -  number  of  bifilar  arms  (N  >  2) 
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Vertical  Linear  Bifilar  BF2 
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The  multi-blade  coordinate  bifilar  arm  flap  displacements 
are  related  to  the  bifilar  arm  physical  flap  displacements  by 

0  =  0O/N  +  <0c  cosi|»  +  0s  sin<J>)(2/N)  (B35) 

See  Figure  9  for  the  positive  directions  of  0  and  connection  node 
displacements.  Figures  B5  to  B7  show  the  mass,  damping,  and 
stiffness  matrices. 


The  applied  load  vector  is 
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The  reaction  load  vector  is 
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( B37  ) 


Figure  2  illustrates  the  positive:  directions  of  the  reaction 
loads . 
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Figures  B8  and  B9  illustrate  the  non-zero  elements  of  the 
mass,  damping,  and  stiffness  matrices  for  the  isolator.  The 
applied  force  vector  is 


{f}  =  {0}  ( B39 ) 

The  reaction  load  vector  is 
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Directions  of  force  components  at  a  connection  node  are  shown  in 
Figure  11. 


Figure  11  illustrates  the  isolator  with  a  single  antiresonant  bar 
weight,  and  Figure  12  is  a  schematic  showing  the  isolator  with  the 
second  antiresonant  bar  weight. 


Properties  of  the  isolator  are  A,  B,  C,  D  -  bar  dimensions  denoted 
by  A^,  Bl,  Cl,  and  D^,  respectively,  in  Figure  12. 
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effective  mass  of  first  isolator  bar  for  yT  dis¬ 
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186 


IP 


lz 


m 


2x 
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w  ,  w2  -  antiresonant  bar  weights 
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ratios  of  natural  frequencies  of  antiresonant 
bars  to  forcing  frequencies  for  yT,  zT,  and  xfi 
displacements,  respectively. 


£ly,  Clz/  C2x  “  damPin9  ratios  of  antiresonant  bars  for  yT,  zT, 
and  x_  displacements,  respectively. 

D 
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“!z>'+<2tlZ  »lz>*l  <B46) 
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I 


zY 


effective  mass  moment  of  inertia  of  second  anti¬ 
resonant  bar  about  y  axis. 


I2Yy/^1-  '“Iy  ^  2+ ^  2^2YUi2Y  ^  2  ^ 


( B47 ) 


1 


1Y' 


mass  moments  of  inertia  of  isolator  bars. 


Stiffnesses  listed  below  are  illustrated  in  Figure  11  where  they 
are  prefixed  with  T. 


K_v,  K„v,  K  equivalent  beam  stif fnessesses  for  anti- 

hi  hZ 

resonant  bar  displacements  x  -  x_,  y  -  y_, 

D  1  O  1 

zfi-zT,  respectively. 

K  K  K 

TX'  TY'  TZ  -  stiffnesses  of  springs  between  bar  pivot  1 

and  connection  node  1. 


KAX'  kay'  kaz 


stiffnesses  of  springs  between  bar  pivot  2 
and  connection  node  2 . 


K0B'  K9T'  K<J>X 


equivalent  beam  stiffnesses  for  antiresonant 

bar  rotations  6„,  0_,  and  , 

B  T  T,  respectively. 


0B  = 

<XT  - 

XB)/C 

( B48 ) 

6T  = 

(ZB  - 

zt)/A 

( B49 ) 

0T  = 

(yB  - 

yt)/a 

( B50 ) 
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Damper  elements  are  in  parallel  with  stiffness  elements  shown  in 
Figures  11  and  12.  Damper  elements  corresponding  to  stiffness 
elements  listed  above  are  obtained  by  replacing  K  by  C.  Damper 
values  are 


where  £ 


TX' 


C, 

c 

c, 

c 

c 

c 

c 

c 

c 

etc . 


TX 

2$tx 

ktx/ujtx 

(  B51  ) 

TY 

= 

2t=Ti 

KTY/uiTY 

( B52  ) 

TZ 

2£tz 

ktz/u)tz 

( B53  ) 

AX 

= 

2^ax 

kax/wax 

( B54  ) 

AY 

= 

2^aykay/ujay 

( B55  ) 

AZ 

= 

2^az 

kaz/u)az 

(B56) 

0B 

= 

2^0B 

K0B/u,0B 

(B57) 

0T 

2^0T 

K0T/w0T 

(B58) 

0T 

= 

2^0B 

K0B/u,0B 

( B59  ) 

• 

are 

damping  ratios,  and 

u)TX,  etc.  are 

frequency  constants  employed  to  define  damper  values. 


The  isolator  with  two  antiresonant  bar  weights  (w^  and  w2 )  is 
referred  to  as  the  3-D  isolator,  and  the  isolator  with  one  anti¬ 
resonant  bar  weight  (w.^)  is  referred  to  as  the  2-D  isolator. 
Properties  of  the  2-D  isolator  are  obtained  from  the  3-D  isolator 
properties  by  equating  to  zero  the  following  parameters. 


W2 '  X2Y'  K0B'  K6T'  C0B'  C0T’ 


The  following  parameters  are  equated  to  one: 


C(=  CL),  D(=  DL). 
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Aeroelastic  Rotor  Type  RE3 


See  Assembly  of  Rotor  Impedance  Matrix  for  coordinates  and 
equations  of  motion. 
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